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XIV 


1.  INTRODUCTION 


The  problem  of  rolled  homogeneous  armor  (RHA)  plate  perforation  by  a 
shaped  charge  jet  is  one  of  ongoing  interest  to  the  U.S.  Army  Ballistic  Research 
Laboratory.  With  the  advent  of  supercomputers,  the  problem  can  now  be  analyzed 
with  more  detail  and  greater  realism  than  has  previously  been  possible.  The  present 
study  presents  an  analysis  of  a  single  case,  that  of  a  particular  shaped  charge 
warhead  fired  at  normal  incidence  and  long  standoff  into  a  13.0-mm-thick  plate. 

The  range  firing  is  described  in  Section  2.  Section  3  follows  with  the  evaluation 
from  experiment  of  three  quantities  pertaining  to  the  perforated  target  plate  that 
provide  convenient  benchmarks  for  hydrocode  calculations.  As  preparation  of  input 
to  the  hydrocode,  the  geometry  of  the  jet’s  leading  particle  is  analyzed  from 
radiographs  in  Section  4.  Features  of  EPIC-2  (Johnson  and  Stryk  1986),  the 
hydrocode  used  in  this  study,  are  reviewed  in  Section  5.  Of  particular  interest  is  its 
eroding  slideline  feature.  Also,  a  series  of  problems  using  only  features  contained 
in  the  1986  version  of  EPIC-2  are  presented  as  an  initial  modeling  effort  in  this 
section.  Section  6  presents  evidence  for  the  failure  mechanisms  of  tensile  voids  and 
shear  bands,  as  observed  by  applying  microscopy  to  specimens  removed  from  the 
perforated  plate.  Section  7  describes  a  model  for  failure  due  to  tensile  voids  that 
has  been  installed  into  EPIC-2.  A  series  of  problems  employing  this  added  feature  is 
discussed.  Section  8  describes  a  model  for  failure  by  shear  bands  which  has  also 
been  installed.  Six  series  of  problems  are  run  using  this  model  in  conjunction  with 
that  for  tensile  voids.  Results  from  these  series  are  compared  with  experiment.  In 
Section  9  one  problem  from  each  of  Sections  5,  7  and  8  is  rerun  with  a  finer  mesh  in 
order  to  observe  the  sensitivity  of  results  to  the  gridding.  Section  10  summarizes 
the  findings  and  suggests  follow-up  work. 

2.  RANGE  FIRING 

The  shaped  charge  warhead  used  has  a  conical  copper  liner  with  a  42°  apex 
angle,  an  81.28-mm  outer-diameter  base,  and  a  1.91-mm  wall  thickness.  The 
explosive  fill  used  is  Composition-B.  It  is  fired  at  a  standoff  of  15.23  charge 
diameters  (C.D.)  into  a  13.0-mm-thick  RHA  plate.  RHA  is  a  medium-carbon, 
quenched  and  tempered,  martensitic  steel.  Its  allowable  ranges  of  chemical 
composition  and  heat  treatment  are  specified  in  U.S.  Department  of  Defense  (1984). 
The  plate’s  cross-section  is  approximately  square,  with  a  197-mm  edge  length.  A 
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Brinell  hardness  number  (BHN)  of  364  for  the  entrance  and  exit  surfaces  is 
obtained.  Figure  1  shows  the  experimental  setup,  including  the  position  and 
orientation  of  four  flash  X-ray  tubes.  The  round  number  is  10771  of  Range  7A. 
Figure  2  displays  two  flash  radiographs  on  one  film  sheet,  taken  at  161.5  and  181.4 
|is  after  explosive  initiation.  Figure  3  shows  the  images  on  the  second  sheet.  These 
apply  at  181.0  and  226.3  jas  after  initiation.  Figure  4  contains  an  enlargement  of  the 
image  at  226.3  ps  after  initiation.  This  time  corresponds  to  about  36  ps  after  intial 
impact  with  the  plate.  Figures  2  through  4  show  the  jet’s  leading  particle  to  have 
completely  perforated  the  plate.  The  trailing  particles  are  aligned  with  the  leading 
particle  and  pass  through  the  hole  that  it  has  made. 

3.  THREE  BENCHMARK  QUANTITIES  EVALUATED  FROM 
EXPERIMENT 

The  perforated  plate’s  entrance  surface  is  shown  in  Figure  5.  Its  exit  surface  is 
presented  in  Figure  6.  The  plate’s  hole  is  reasonably  circular,  except  for  a  single 
notched  region.  Away  from  this  notch,  the  hole  is  assumed  to  have  been  created 
entirely  by  the  jet’s  leading  particle.  This  assumption  is  supported  by  the  degree  of 
alignment  in  the  jet’s  particles  that  was  noted  in  the  radiographs  above.  The  hole  is 
fitted  by  a  circle  (Figure  7),  and  radii  spaced  at  25.7°  are  constructed  (Figure  8). 
Slices  of  plate  material  corresponding  to  five  adjacent  radii  are  removed  (Figures  9 
and  10).  The  boundary  of  the  uppermost  slice  in  Figure  10  is  digitized  to  describe 
the  hole’s  profile  (Figure  11). 

r^fe  is  defined  as  the  value  determined  from  experiment  of  the  final  hole 
boundary’s  radius  averaged  over  the  plate’s  thickness.  The  average  radius  of  the 
red  surface  in  Figure  11  is  determined  by  numerical  integration  of  the  digitized 
coordinates.  The  result  is  a  value  of  18.2  mm  for  r^xfe ,  so 

fhofe  =  18.2  mm  (1) 

This  quantity  will  serve  as  a  convenient  benchmark  by  which  to  evaluate  hydrocode 
solutions.  Two  other  quantities  that  will  also  serve  this  function  are  AMexp '  and 
t eg$' .  A Mexp  is  defined  to  be  the  experimentally-determined  total  mass  lost  by  the 
target  plate.  tegf  is  the  experimentally-determined  time  after  initial  impact  at  which 
the  plate  has  lost  95%  of  this  mass.  No  data  by  which  to  evaluate  A Mexp  and  t 
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were  obtained  from  Round  (Rd.)  10771. 


Fortunately,  two  more  recent  experiments,  Rds.  4189  and  4190  of  Range  16,  do 
provide  such  data.  In  both  of  these  experiments,  the  same  shaped  charge  warhead 
employed  in  Rd.  10771  was  fired  into  an  RHA  plate  at  a  standoff  of  975  mm  (12.0 
C.D.)  The  impact  conditions  at  this  standoff  are  quite  similar  to  those  at  15.2  C.D. 
In  both  cases  the  entire  jet  has  broken  up  into  particles  that  have  probably  ceased  to 
stretch.  Spacings  between  successive  particles  are  greater  at  15.2  C.D.  than  at  12.0 
C.D.,  but  Figures  3  and  4  establish  that  the  fastest  particles  and  probably  most 
particles  pass  through  the  hole  made  by  the  first  particle  without  contacting  the 
plate.  The  target  plates  in  Rds.  4189  and  4190  had  respective  thicknesses  of  12.6 
and  12.8  mm,  both  close  to  the  13.0  mm  value  of  Rd.  10771.  Thus  data  from  Rds. 
4189  and  4190  are  assumed  transferrable  to  Rd.  10771. 

In  Rd.  4190  the  plate  was  weighed  before  and  after  the  firing  to  produce  a 
measurement  of  136  g  for  A Mexp  ,  so 

A MexP  =  136  g  (2) 

A  total  of  six  flash  radiographs  was  obtained  from  Rds.  4190  and  4189.  As  shown 
in  Figure  12,  these  radiographs  visualize  the  debris  patterns  at  1.7,  18.2,  33.2,  73.9, 
143.2  and  213.0  ps  after  initial  impact.  The  three  earliest  views  were  obtained  from 
Rd.  4190,  and  the  three  latest  from  Rd.  4189.  The  first  RHA  fragments  to  form  are 
small  and  are  seen  to  keep  pace  with  the  fastest  moving  jet  particles.  Much  larger 
fragments  appear  at  later  times.  In  Section  8  distinct  micromechanical  mechanisms 
will  be  suggested  for  these  two  classes  of  fragments.  It  seems  from  Figure  12  that 
substantial  material  has  separated  from  the  plate  subsequent  to  73.9  ps  after  impact. 
A  rough  estimate  for  t e$  is  therefore 

te9f  >  73.9  ps  (3) 

The  difficulty  with  assigning  a  time  of  fragment  separation  from  the  target  plate  on 
the  basis  of  radiographs  such  as  Figure  12  should  be  borne  in  mind,  however.  Such 
radiographs  directly  show  only  the  later  time  at  which  the  separated  fragment  has 
moved  a  noticeable  distance  from  the  plate. 
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4.  LEADING  JET  PARTICLE  DESCRIPTION 


Figure  2  shows  the  leading  jet  particle  at  161.5  and  at  181.4  ps  after  initiation. 
These  times  correspond  to  28.5  and  8.6  ps  before  impact,  respectively.  Based  on 
these  two  images,  the  leading  particle’s  impact  velocity  is  determined  to  be  7.73 
mm/ps. 

A  cylindrical  coordinate  system  is  introduced,  with  its  origin  at  the  center  of  the 
target  plate’s  entrance  surface.  Axial  coordinate  z  measures  distance  from  this 
surface,  with  positive  values  in  the  region  between  the  plate  and  the  warhead’s 
original  position.  The  digitized  image  of  the  jet’s  leading  particle  at  8.6  ps  before 
impact  is  mapped  onto  this  cylindrical  coordinate  system  by  translation  and  rotation. 
The  mapping  is  defined  by  three  constraints:  (i)  and  (ii)  The  particle’s  volume  and 
centroid  location,  when  computed  based  on  the  portion  on  one  side  of  the  r=  0 
centerline,  agree  with  those  computed  using  the  portion  on  the  other  side  of  the 
centerline,  (iii)  Following  translation  and  rotation,  the  point  on  the  digitized 
boundary  that  intersects  the  centerline  at  the  lowest  z  coordinate  is  assigned  the 
position  z=  0,  i.e.,  is  placed  in  contact  with  the  center  of  the  plate’s  entrance 
surface.  In  addition  to  producing  a  geometry  on  which  to  base  input  to  an 
axisymmetric  hydrocode  analysis,  the  mapping  determines  the  particle’s  length, 
volume  and  centroid  position  to  be  20.6  mm,  353.2  mm3,  and  7.97  mm  above  the 
base,  respectively. 

A  five-parameter  model  for  the  leading  particle’s  geometry  is  shown  in  Figure 
13.  The  model’s  length  Lp,  volume  Vp,  and  centroid  position  zp  are  expressed  in 
terms  of  its  five  parameters  by 


Lp  =  Lc  +  Ls  +  Lg 


(4) 


Vp  -  nRs* 


1 


Lc  +  Ls  +  LgX 


(5) 
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(6) 


^  Lc  LB  +  Ls  +  ^Lc  +  Ls  LB  +  ^  Ls  +  2  ^b2X 


~LC  +  Ls  +  Lbx 


where 


1  2_  &b  ,  ,  2.  Rb_  1 

y  =  1  +  ~  1  +  —  —  -  1 

2  Rs  3  I** 


LP,  Vp  and  zp  are  assigned  the  values  20.6  mm,  353.2  mm3,  and  7.97  mm, 
respectively,  as  evaluated  above  on  the  basis  of  a  radiograph.  Equations  (4),  (5) 
and  (6)  can  then  be  viewed  as  three  equations  in  the  five  unknowns  Lc,  Ls,  LB,  Rs 
and  Rb.  A  solution  is  obtained  for  Lc,  Ls  and  LB  in  terms  of  Rs  and  RB. 
Corresponding  to  a  guess  for  the  values  of  Rs  and  RB,  the  quantities  Lc,  Ls  and  LB 
are  evaluated  and  the  model  of  Figure  13  is  drawn.  This  procedure  is  applied 
iteratively  until  a  close  visual  fit  to  the  jet  particle’s  shape  according  to  the 
radiograph  is  obtained.  The  image  in  Figure  14,  which  corresponds  to  the  values 
for  Lc,  Ls,  Lb,  Rs  and  RB  that  are  listed  in  Table  1,  is  deemed  to  be  an  acceptable  fit 
to  the  radiograph. 

5.  EPIC-2  SIMULATIONS  WITHOUT  FAILURE  MODELING 

5.1  Features  of  EPIC-2.  The  1986  version  of  Lagrangian  hydrocode  EPIC-2 
(Johnson  and  Stryk  1986)  is  used  to  model  the  Rd.  10771  experiment.  The  two 
metals,  copper  and  steel,  are  both  handled  in  the  same  fashion.  Their  dilatational 
deformation  is  governed  by  the  Mie-  Grimeisen  equation  of  state 

'  - 

p=  (*,*!+  K2p2+  K3[i3]  1-TTp  +  r(l  +  n)e  (8) 


where 


\i  =  -  1 

Po 
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p  is  pressure,  or  negative  hydrostatic  stress,  p  is  current  density,  po  is  undeformed 
density,  e  is  the  internal  energy  per  undeformed  volume,  and  T  is  the  Gruneisen 
coefficient.  Values  for  Kx,  K2>  K3  and  T  are  presented  for  OFHC  copper  and  for 
304  stainless  steel  in  Kohn’s  handbook  (Kohn  1965).  These  parameter  values  are 
displayed  in  Table  2,  in  which  Kohn’s  values  for  austenitic  304  stainless  steel  have 
been  applied  to  martensitic  RHA. 


The  distortional  deformation  of  copper  and  steel  is  treated  by  means  of  a 
plasticity  model,  which  applies  the  von  Mises  yield  condition  as  follows:  At  each 
cycle,  deviatoric  Cauchy  stress  components  srr,  s sqq  and  srz  are  computed  by 
applying  elastic  shear  modulus  G  to  the  total  strain  increment.  Equivalent  stress  a 
is  then  computed  using 


<5  = 


'zz 


+  see  )  +  3s 


2 


n 


(10) 


and  is  compared  to  the  current  flow  stress  ay  ,  a  measure  of  the  state  of  the 
material.  If  a  <  ay  ,  only  elastic  deformation  occurs  during  that  cycle  and  the 
computed  stress  components  are  unaltered.  If  o  >  ay ,  plastic,  non-recoverable 
deformation  occurs  as  well  as  elastic.  If  a>  ay  ,  stresses  are  adjusted  along  a  radial 
path  in  stress  space  so  that  o  becomes  equal  to  oy . 


Flow  stress  oy  is  computed  using  the  Johnson-Cook  strength  model  (Johnson 
and  Cook  1983).  According  to  this  model,  ay  is  a  function  of  equivalent  plastic 
strain  ep ,  equivalent  plastic  strain  rate  ep,  and  homologous  temperature  T*.  The 
functional  form  is 


ov  = 


1  +  C  In 


1-  ( T*)m 


(11) 


where 


ip  = 


~DP  Dp 
2  uiju  ij 


(12) 


and 


(13) 


zp  =  \  kpdt 
o 


t  is  time  and  Df;  is  the  plastic  component  of  the  rate  of  deformation  tensor.  T*  is 
related  to  room  temperature  Tr  and  melting  temperature  Tm  by 


T * 


T-  Tr 
Tm~  Tr 


(14) 


where  T  is  the  current  temperature,  computed  using  a  thermodynamic  relation. 
Equation  (11)  expresses  strain  hardening,  strain  rate  hardening,  and  thermal 
softening  as  separate  factors.  Hardening  is  assumed  to  be  isotropic.  Constants  G, 
A,  B,  N,  C,  M,  7V,  and  Tm  are  assigned  their  values  presented  for  OFHC  copper  and 
for  4340  steel  by  Johnson  and  Cook  (1983).  These  values  are  listed  in  Table  2. 

The  slideline  erosion  algorithm  of  EPIC -2  enables  the  application  of  this 
Lagrangian  hydrocode  to  problems  in  penetration  mechanics.  A  description  of  this 
essential  feature  is  reserved  for  the  next  subsection,  5.2. 


In  addition  the  1986  version  of  EPIC-2  contains  a  failure  model,  or  algorithm 
to  cause  stress  reduction  within  a  finite  element  subsequent  to  the  local  satisfaction 
of  some  onset  criterion.  This  Johnson-Cook  failure  model  is  phenomenological  and 
was  derived  >  fit  data  from  Hopkinson  bar  tests  and  quasi-static  tensile  tests 
(Johnson  and  Cook  1985).  Such  tests  involve  relatively  low  strain  rates,  so  the 
model  is  judged  to  be  inappropriate  for  the  present  application.  In  particular 
descriptions  in  Johnson  and  Cook  (1985)  of  damaged  specimens  contain  no  mention 
of  shear  banding,  which  in  Section  6  will  be  shown  to  constitute  one  of  two 
micromechanical  mechanisms  for  material  failure  found  to  be  operative  in  RHA 
specimens  from  Rd.  10771.  The  other  micromechanical  mechanism  is  ductile  voids, 
which  apparently  was  also  active  in  the  experiments  of  Johnson  and  Cook. 
Nevertheless,  the  range  of  ambient  conditions  on  strain  rate  and  temperature  were 
indeed  much  smaller  in  their  experiments  than  in  Rd.  10771,  so  a  much  simpler 
model  for  void  failure  will  be  utilized  instead  in  this  report. 

5.2  Problem  Series  I.  The  axisymmetric  mesh  used  in  all  calculations  of  this 
section  as  well  as  Sections  7  and  8  is  shown  in  Figure  15.  It  consists  of  three-node 
triangular  elements.  The  copper  projectile’s  mesh  closely  approximates  the 
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geometry  of  the  five-parameter  model  with  parameter  values  listed  in  Table  1.  It 
contains  124  node  points  and  202  elements.  The  target’s  mesh  approximates  the 
square  steel  plate  by  a  circular  disk  with  a  thickness  of  13.0  mm  and  a  radius  of 
104.0  mm.  Individual  elements  in  the  target  are  isosceles  right  triangles,  the  longest 
edge  of  which  has  a  length  of  1.3  mm.  This  corresponds  to  ten  layers  of  crossed 
triangles  through  the  target  plate’s  thickness.  The  target’s  mesh  includes  3200  of 
these  elements  and  1691  node  points. 

The  mesh  contains  two  slidelines,  both  coincident  with  the  projectile/target 
interface.  One  slideline  has  the  projectile’s  surface  nodes  serving  as  master  nodes 
and  the  target’s  surface  nodes  as  slaves.  In  the  other  slideline  this  assignment  is 
reversed. 

The  erosion  feature  is  operative  for  both  slidelines.  This  allows  an  element 
having  one  or  more  corner  node  on  the  master  surface  of  the  slideline  to  be 
discarded  from  the  problem  in  the  sense  that  all  stresses  in  the  element  are 
thereafter  set  to  zero.  Since  both  copper  projectile  nodes  and  steel  target  nodes 
alternately  serve  as  master  nodes  for  one  of  the  two  slidelines,  the  procedure 
provides  a  means  for  modeling  both  projectile  erosion  and  target  hole  formation. 
EPIC-2  uses  three  criteria  to  trigger  erosion  of  a  given  surface  element.  One 
criterion  is  a  user-supplied  cutoff  value  for  equivalent  plastic  strain,  strode-  Another 
is  an  angle  cutoff,  whereby  none  of  the  element’s  vertex  angles  is  allowed  to  become 
less  than  a  certain  acute  value.  These  two  criteria  limit  the  degree  of  distortion  in 
the  element.  The  third  criterion  discards  a  group  of  one  or  more  elements  that  is 
connected  to  the  master  surface  by  a  single  vertex  node.  This  criterion  is  motivated 
by  the  need  to  prevent  such  elements  from  crossing  the  slideline  surface. 

Following  erosion  of  a  slideline  element,  its  three  associated  node  points,  to 
which  its  entire  mass  has  been  lumped,  remain  in  the  problem.  One  or  two  may 
become  effectively  disembodied  from  the  remaining  projectile  or  target.  When  this 
occurs,  those  disembodied  nodes  are  converted  from  master  into  slave  nodes.  They 
subsequently  interact  with  the  remaining  master  surface  and  participate  in 
momentum  transfer  across  the  slideline.  Furthermore,  previously  interior  nodes  that 
now  find  themselves  located  on  the  master  surface  become  identified  as  the 
slideline’s  master  nodes.  EPIC-2’s  eroding  slideline  algorithm  is  described  further 
in  Stecher  and  Johnson  (1984). 
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In  later  sections  of  this  report,  two  failure  models  will  be  developed  to  cause 
stress  reduction  in  el  uents  not  necessarily  located  along  the  current  slideline. 

Before  considering  this  additional  feature,  consider  a  series  of  five  problems  that  are 
run  using  only  those  features  already  discussed.  It  will  be  referred  to  as  Problem 
Series  I.  Erosion  strain  zperode  is  assigned  the  values  0.25,  0.50,  0.75,  1.25  and  2.00, 
as  listed  in  Table  3.  Each  problem  begins  at  the  instant  of  initial  impact  and  is  run 
until  500  ps  after  impact.  Mesh  plots  for  Problem  I-C,  for  which  e^rode  equals  0.75, 
are  shown  at  3.01,  4.02,  and  500.04  jj. s  after  impact  in  Figures  16,  17,  and  18, 
respectively.  Perforation  has  occurred  between  3.01  and  4.02  ps  after  impact.  This 
implies  a  penetration  velocity  between  3.23  and  4.32  mm/ps  averaged  over  the 
course  of  perforation. 

Figure  19  shows  the  mesh  from  Problem  I-C  at  10  |is  after  impact,  with  the 
nodes  disembodied  by  slideline  erosion  indicated.  Each  of  these  nodes  is  assigned  by 
means  of  its  color  to  one  of  the  four  categories:  "target  debris",  "projectile  debris”, 
"target  ejecta",  and  "projectile  ejecta".  "Target  debris"  denotes  disembodied  nodes 
of  RHA  material  having  a  negative  axial  velocity  z.  According  to  the  coordinate 
system  of  Figure  19,  nodes  with  negative  axial  velocity  travel  in  the  same  direction 
as  the  jet  particle  and  away  from  the  original  (before  detonation)  warhead  position. 
"Target  ejecta"  denotes  disembodied  nodes  of  RHA  material  having  a  positive  axial 
velocity,  or  traveling  towards  the  original  warhead  position.  Similarly,  "projectile 
debris"  and  "projectile  ejecta"  denote  disembodied  nodes  of  copper  material  with  a 
negative  and  a  positive  axial  velocity,  respectively. 

Now  consider  results  for  the  average  hole  radius  at  500  ps  after  initial  impact, 
rhoie .  Figure  20  redraws  the  boundary  of  the  target  mesh  of  Figure  18.  The  portion 
of  this  boundary  associated  with  the  perforation  hole  boundary  is  shown  in  red. 

^hoie  is  computed  from  Figure  20  by  numerical  integration  of  nodal  coordinates 
along  this  red  section.  The  result  for  Problem  I-C  is  11.0  mm,  substantially  smaller 
than  the  experimental  result  of  18.2  mm  obtained  from  Figure  11.  Results  for  rhote 
from  Problem  I-C  and  the  other  four  problems  in  the  series  are  listed  in  Table  3  and 
plotted  vs.  E^rode  in  Figure  21.  The  computed  values  for  rhole  are  seen  to  increase 
with  increasing  zperode ■  However,  even  for  an  zperode  of  2.00  (Problem  I-E),  the 
computed  rhoie  of  13.9  mm  is  still  substantially  smaller  than  the  experimental  result. 
The  mesh  of  the  target  plate  at  500  ps  for  Problem  I-E  is  presented  in  Figure  22. 
The  large  value  of  zperode  f°r  this  problem  has  allowed  highly  distorted  elements  to 
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be  retained  in  the  final  mesh.  This  consideration  leads  to  the  decision  not  to  apply 
values  of  £perode  in  excess  of  2.00,  even  though  the  trend  in  Figure  21  suggests  that 
such  values  might  yield  more  realistic  results  for  rhoie. 

Figure  23  plots  the  total  mass  that  has  eroded  from  the  target  plate  vs.  time 
after  impact  for  the  five  problems  of  Problem  Series  I.  These  total  masses  are 
computed  for  a  given  time  step  by  summing  the  masses  of  all  steel  elements  that 
have  eroded  at  or  before  that  time  step.  Results  for  AM,  the  total  mass  lost  at  500 
ps,  and  /  95,  the  time  at  which  the  plate  has  lost  0.95  AM  of  mass,  can  be  obtained 
from  Figure  23.  These  results  are  collected  in  Table  3.  AM  results  are  rather 
insensitive  to  £%rode  for  zperode  equal  to  0.25,  0.50,  and  0.75.  However,  by  an  £%rotje 
of  1.25,  AM  is  seen  to  decrease  with  increasing  £perode ■  This  last  trend  is  in  contrast 
to  the  increase  in  hole  radius  rhoie  with  increasing  zperPje  that  was  noted  above.  This 
discrepancy  between  relationships  of  rhole  and  AM  to  is  resolved  by  a 
comparison  of  the  mesh  plots  in  Figures  18  and  22.  Hole  growth  is  seen  to  have 
been  accompanied  by  the  removal  of  fewer  RHA  elements  in  Problem  I-E  than  in 
Problem  I-C.  Instead,  in  Problem  I-E  many  highly  distorted  elements  have  been 
retained,  leading  to  much  radial  compression  of  material  near  the  hole  "throat",  or 
most  narrow  cross-section,  and  to  the  formation  of  prominent  lips  at  the  entrance 
and  exit  surfaces  along  the  hole  boundary.  AM  values  for  the  five  problems  are  all 
much  smaller  than  the  experimental  value  of  136  g  in  Equation  (2).  This  is 
consistent  with  the  fact  that  computed  rhoie  values  are  all  smaller  than  the 
measurement  from  Rd.  10771.  Results  for  t  95  in  Table  3  are  all  substantially 
smaller  than  the  experimental  lower  bound  of  Equation  (3). 

The  results  from  Problem  Series  I  can  be  summarized  as  follows:  For  the  mesh 
of  Figure  15  and  for  the  parameter  £per0de  varied  simultaneously  for  both  slidelines  in 
the  range  of  0.25  to  2.00,  (i)  rhoie  is  smaller  than  the  measurement  from  Rd.  10771, 
(ii)  AM  is  smaller  than  the  measurement  from  Rd.  4190,  and  (iii)  t  95  is  smaller  than 
the  range  indicated  by  Rds.  4190  and  4189.  In  short  the  computed  holes  are  too 
small,  and  their  growths  are  completed  prematurely. 
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6.  EXPERIMENTAL  OBSERVATIONS  OF  VOIDS  AND  SHEAR  BANDS 


At  this  point  the  perforated  target  plate  of  Rd.  10771  is  studied  further  to  assess 
the  micromechanical  processes  that  it  has  undergone.  RHA  specimens  are  extracted 
from  the  perforated  plate  in  the  vicinity  of  the  hole.  They  are  then  polished,  etched 
in  2%  Nital,  and  examined  using  photomicroscopy  and  scanning  electron 
microscopy.  Figure  24  displays  a  fracture  surface  that  has  been  formed  by  void 
coalescence.  Figure  25  shows  a  shear  band  with  a  characteristic  6-pm  width,  while 
Figure  26  displays  a  crack  that  runs  along  part  of  a  shear  band’s  length.  From  these 
figures,  it  is  concluded  that  fracture,  and  therefore  fragmentation  of  RHA,  are 
associated  with  the  two  phenomena  of  voids  and  shear  bands  that  result  from 
interaction  with  the  jet. 

The  extracted  slice  of  RHA  target  plate  from  Rd.  10771  that  was  photographed 
to  produce  Figure  11  is  now  polished  and  examined  with  an  optical  microscope  at 
magnifications  up  to  1000X.  The  goal  is  to  observe  the  locations  of  prominent  voids 
and  shear  bands  in  the  final  target  plate  cross-section.  Here  "prominent"  is  used  to 
mean  "observable  with  an  optical  microscope  at  1000X  magnification."  The  results 
are  sketched  in  Figure  27.  Voids  and  shear  bands  are  seen  to  be  concentrated  in  the 
hole  boundary’s  region  that  is  located  near  the  plate’s  midsurface  and  that  protrudes 
into  the  hole  to  form  the  throat.  At  least  nine  large,  discrete  voids  and  eleven  shear 
bands  can  be  observed  in  this  region.  Each  shear  band  contacts  either  a  void,  the 
hole  boundary,  or  in  some  cases  both.  Two  or  three  prominent  shear  bands  are  also 
observed  near  the  intersection  of  the  hole  boundary  with  the  plate’s  entrance 
surface. 

7.  EPIC-2  SIMULATIONS  WITH  THE  TENSILE  FAILURE  MODEL 

7.1  The  Tensile  Failure  Model.  The  voids  in  the  RHA  plate  noted  in  Figure  24 
are  assumed  to  have  nucleated  and  grown  under  the  action  of  a  tensile  stress  field. 

A  tensile  failure  model  to  simulate  effects  of  these  voids  at  the  finite  element  level  is 
now  developed  and  added  to  EPIC-2.  In  contrast  to  the  slideline  erosion  algorithm, 
this  tensile  failure  model  will  be  applied  to  target  and  jet-particle  elements 
throughout  the  domain,  both  on  and  off  the  slideline.  The  model  consists  of  two 
parts:  (i)  a  failure  onset  criterion  and  (ii)  a  prescription  for  post-onset  stress 
reduction. 
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The  onset  criterion  of  the  tensile  failure  model  is  developed  first.  The 
approximately  spherical  shape  of  individual  voids  in  the  photomicrograph  of  Figure 
24  suggests  that  their  growth  was  driven  by  the  spherically  symmetric  portion  of  the 
stress  tensor,  namely  the  hydrostatic  stress.  Pressure  is  negative  hydrostatic  stress, 
defined  in  terms  of  arr,  ozz  and  Oqq,  the  radial,  axial  and  circumferential  normal 
components  of  the  Cauchy  stress  by 


P 


1 

3 


+  <*zz 


+ 


(15) 


The  RHA  specimen  of  Figure  24  was  removed  from  the  remaining  target  plate 
in  the  vicinity  of  the  perforation  hole.  The  figure  therefore  displays  material  at 
about  19  mm  from  the  shot  line.  The  observation  just  made  regarding 
spheroidally-shaped  voids  is  assumed  to  apply  throughout  the  RHA  target  plate. 

This  last  assumption  is  challengeable  in  that  during  the  perforation  process,  RHA 
material  that  was  originally  located  closer  to  the  shot  line  presumably  experienced 
much  larger  values  for  such  field  variables  as  strain,  strain  rate  and  temperature 
than  did  the  material  of  Figure  24.  However,  the  assumption  receives  some  support 
from  publications  by  researchers  at  SRI  International  (e.g.,  Barbee  et  al.  1972, 
Seaman  et  al.  1976).  These  publications  have  established  that  this  general  situation 
of  fractures  formed  by  coalescing  spheroidal  voids  occurs  quite  commonly  in  shock 
loaded  metals.  They  call  such  fracture  "ductile",  as  opposed  to  their  category  of 
"brittle"  fracture,  which  is  characterized  by  planar  cracks.  In  particular  Barbee  et  al. 
(1972)  has  shown  the  situation  to  be  at  least  sometimes  applicable  to  OFHC  copper. 
On  this  basis  a  pressure  criterion  for  tensile  failure  onset  will  be  applied  to  the 
copper  jet  particle  as  well  as  to  the  RHA. 

In  an  element  of  either  RHA  target  material  or  OFHC  copper  jet  material,  the 
onset  of  tensile  failure  is  assumed  to  occur  when  the  pressure,  which  at  a  given  time 
step  is  constant  in  an  element,  becomes  less  than  or  equal  to  a  negative-valued 
material  parameter  Pfan,  i.e.. 


p  <,  pfail  <0  (16) 

For  each  metal,  pfaU  is  taken  to  be  a  constant,  independent  of  both  the  local  current 
values  and  the  histories  of  all  field  variables. 
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The  prescription  for  post-onset  behavior  will  be  considered  next.  Once  the 
failure  onset  criterion  of  Equation  (16)  is  met  in  an  element,  that  element  is  no 
longer  allowed  to  support  hydrostatic  tensile  stress  or  deviatoric  stresses.  The 
element  is  only  allowed  to  support  hydrostatic  compression.  This  involves 
modifications  to  the  EPIC-2  hydrocode  which  at  each  time  step  affect  all  elements 
that  have  either  then  or  at  some  previous  time  step  satisfied  Equation  (16).  For  all 
such  elements,  deviatoric  stresses  are  not  computed,  and  positive  hydrostatic  stress 
values  that  result  from  evaluation  of  the  Mie-  Griineison  equation  of  state  are 
altered  and  set  to  zero. 

Note  that  this  post-onset  prescription  involves  the  instantaneous  reduction  to 
zero  of  certain  stress  components  once  the  tensile  failure  onset  criterion  has  been 
met.  The  prescription  therefore  ignores  the  possiblity  for  gradual  degradation  of  a 
material’s  load-bearing  capability  prior  to  its  complete  local  failure.  This  neglect 
may  be  justifiable  in  problems  involving  shock  loading,  since  the  time  scale  of 
material  degradation  may  be  much  smaller  than  those  which  pertain  to  other  aspects 
of  the  problem. 

The  complete  tensile  failure  model,  including  the  onset  criterion  and  the  post¬ 
onset  prescription,  introduces  a  single  additional  input  parameter  for  each  material, 
namely  pfan.  The  issue  to  be  addressed  next  is  the  assignment  of  a  value  to  pfaii  for 
each  of  the  two  materials.  Relevant  data  were  obtained  by  Rinehart  (1951)  from 
spallation  experiments  performed  on  several  metals,  including  copper  and  4130 
steel.  Spallation  is  defined  as  material  failure  caused  by  the  tensile  stress  field  that 
results  from  reflection  of  a  compressive  pulse  from  a  free  surface.  In  each 
experiment  a  cylindrical  explosive  charge  was  placed  on  one  surface  of  a  metal 
plate.  A  pellet  composed  of  the  same  material  as  the  plate  was  placed  on  the 
opposite  surface  to  serve  as  a  momentum  trap.  For  each  metal  studied,  Rinehart 
determined  a  value  for  the  "critical  normal  fracture  stress".  This  is  the  value 
attained  in  the  region  of  spall  fracture  by  the  axial  normal  stress  component.  The 
values  that  Rinehart  reported  are  2.96  GPa  for  copper  and  3.03  GPa  for  4130  steel. 
If  one  assumes  that  a  condition  of  uniaxial  strain  applied  to  his  experiments,  the 
pressure  corresponding  to  each  critical  axial  normal  stress  can  be  computed  by 
subtracting  2ay/3  from  the  reported  value,  where  ay  is  the  local  flow  stress  (Zukas 
1982).  The  flow  stress  at  initial  yield  is  determined  in  Johnson  and  Cook  (1983)  to 
be  0.09  GPa  for  OFHC  copper  and  0.79  GPa  for  4340  steel.  If  one  also  assumes 
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that  the  flow  stress  levels  attained  in  Rinehart’s  experiments  were  never  caused  by 
strain  or  strain-rate  hardening  to  greatly  exceed  these  values,  a  reasonable  value  for 
pfaii  is  found  to  be  -3.00  GPa  (-30.0  kbar)  for  both  OFHC  copper  and  RHA  steel. 
This  value  will  be  applied  to  both  materials  in  all  calculations  that  follow. 

The  tensile  failure  model  that  has  been  developed  can  be  summarized  as 
follows:  In  each  element  of  either  RHA  target  or  OFHC-copper  jet-particle 
material,  tensile  failure  first  occurs  when  Equation  (16)  is  satisfied.  At  that  time 
step  and  at  all  time  steps  thereafter,  the  element  is  not  allowed  to  support 
hydrostatic  tensile  stress  or  any  component  of  deviatoric  stress.  Only  hydrostatic 
compression  is  allowed.  The  model  introduces  the  single  additional  material 
constant  pfaih  which  is  assigned  the  value  of  -3.00  GPa  for  both  RHA  and  OFHC 
copper. 

7.2  Problem  Series  II.  A  second  series  of  five  problems.  Problem  Series  II,  is 
run  with  this  tensile  failure  model  operative.  Throughout  the  series,  pfau  is  assigned 
the  value  of  -3.00  GPa  for  both  RHA  and  copper.  As  in  Problem  Series  I,  erosion 
strain  e£rode  is  again  varied  simultaneously  for  both  slidelines  over  the  range  of  0.25 
to  2.00.  Mesh  plots  for  Problem  II-C,  with  zperode  set  to  0.75,  are  shown  at  1.02, 
3.02  and  500.00  ps  after  impact  in  Figures  28,  29  and  30,  respectively.  In  these 
figures  elements  that  are  colored  blue  are  those  that  have  failed  in  tension  according 
to  the  model.  The  phenomenon  of  spallation,  or  tensile  failure  brought  on  by 
reflection  of  a  compressive  shock  wave  at  a  free  surface,  is  exhibited  in  Figure  29. 
Spallation  is  predicted  to  occur  prior  to  perforation  and  in  a  small  region  directly  in 
the  path  of  the  leading  jet  particle.  Figure  30  shows  a  single  non-eroded  element  in 
the  final  target  plate,  located  near  the  plate’s  midsurface,  in  which  tensile  void 
failure  has  occurred.  This  elevation,  or  z  coordinate,  of  void  location  is  in  some 
agreement  with  experiment  in  that  all  voids  in  Figure  27  reside  within  the 
midsurface  protuberance. 

Now  consider  final  mesh  plots  from  the  other  problems  in  this  series.  Figures 
31  through  34  contain  mesh  plots  at  500  ps  from  Problems  II-A,  II-B,  II-D  and 
II-E,  respectively.  No  failed  elements  are  present  in  these  final  mesh  plots  from 
Problems  II-A  and  II-B.  In  Problem  II-D  all  failed  elements  are  located  near  the 
hole  boundary  and  slightly  to  the  exit  side  of  the  midsurface.  A  larger  portion  of 
the  cross-section  has  failed  in  tension  in  Problem  II-D  than  in  Problem  II-C. 
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Moreover,  the  pattern  of  failed  elements  from  Problem  II-D  is  in  better  agreement 
with  the  experimental  distribution  of  Figure  27  than  is  the  pattern  in  Figure  30  from 
Problem  II-C.  Also  note  in  Figure  33  that  in  Problem  II-D  the  higher  value  of  tperode 
has  allowed  highly  distorted  RHA  elements  to  escape  erosion.  The  result  is 
prominent  lips  on  the  entrance  and  exit  surfaces  at  the  hole  boundary.  Such  lips  are 
not  present  on  the  experimental  cross-section  in  Figures  11  and  27.  These  two 
figures  were  obtained  from  the  top  slice  shown  in  Figure  10.  If  instead  one  of  the 
three  middle  slices  in  Figure  10  were  focused  on,  distinct  lips  would  have  been 
included  in  the  experimental  profile. 

In  Problem  II-E  voids  are  present  in  a  still  larger  portion  of  the  cross-section 
than  in  Problem  II-D.  The  lips  around  the  hole  boundary  at  the  entrance  and  exit 
surfaces  have  become  larger  as  well  in  Problem  II-E  and  are  now  clearly  larger  than 
in  any  of  the  five  experimental  slices  in  Figure  10. 

Results  for  rhole,  AM  and  f  95  obtained  from  the  five  problems  of  Problem 
Series  II  are  listed  in  Table  4.  rhoU  is  plotted  as  a  function  of  zperode  for  this  series  in 
Figure  35.  The  mass  lost  from  the  plate  vs.  time  is  plotted  for  the  series  in  Figure 
36. 


A  comparison  of  Table  4  and  Figures  35  and  36  with  Table  3  and  Figures  21 
and  23  shows  that  in  general  little  change  in  computed  rhoie,  AM  and  t  95  values  has 
accompanied  the  addition  of  the  tensile  void  failure  model  with  the  value  of  -3.00 
GPa  assigned  to  Pfaii  for  both  copper  and  RHA.  For  all  five  problems  in  Problem 
Series  II,  rhoie  and  AM  are  smaller  than  the  experimental  values  in  Equations  (1)  and 
(2),  respectively.  Thus  the  computed  holes  are  still  too  small.  For  Problems  II-A, 
II-B,  II-C  and  II-D,  t  95  is  also  much  smaller  than  the  experimental  result  in 
Equation  (3),  indicating  that  hole  formation  has  completed  prematurely.  Problem 
II-E  is  exceptional  in  that  it  does  yield  a  t  95  value  consistent  with  the  experimental 
result.  Sensitivity  of  the  results  to  the  values  used  for  pfau  for  the  two  materials  has 
not  been  studied. 
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8.  EPIC-2  SIMULATIONS  WITH  THE  TENSILE  AND  SHEAR  BAND 


FAILURE  MODELS 

8.1  The  Shear  Band  Failure  Model.  The  algorithm  inserted  into  EPIC-2  to 
model  material  failure  from  shear  bands  contains  the  same  two  generic  parts  as  does 
the  tensile  failure  model.  These  are  a  failure  onset  criterion  and  a  prescription  for 
post-onset  stress  reduction.  The  criterion  adopted  is  called  "Zener-Hollomon",  after 
its  apparent  originators  (Zener  and  Hollomon  1944).  In  a  given  time  step  the 
product  [A+  is  evaluated  for  each  RHA  element  that  has  not  yet 

failed  or  eroded  and  that  has  undergone  further  plastic  flow  during  that  time  step. 
(The  shear  band  model  is  not  applied  to  the  copper.)  This  product  includes 
contributions  to  the  flow  stress  by  plastic  straining  and  by  temperature.  If  the 
element  also  underwent  plastic  flow  during  the  immediately  preceding  time  step  and 
if  the  above  product  has  decreased  from  its  level  during  that  preceding  time  step, 
then  the  element  is  deemed  to  have  begun  to  fail  by  shear  banding.  In  effect,  the 
criterion  compares  the  strength  increase  due  to  plastic  straining  (work  hardening)  to 
its  decrease  due  to  temperature  rise  (thermal  softening).  Failure  occurs  when  the 
latter  exceeds  the  former. 

Once  this  onset  criterion  has  been  satisfied  in  an  element,  at  that  and  all 
subsequent  time  steps,  deviatoric  stresses  are  reduced  from  their  computed  values. 
However,  hydrostatic  tension  is  not  altered,  in  contrast  to  the  post-onset  modeling 
for  the  case  of  tensile  failure.  Another  difference  is  that  the  deviatoric  stresses  are 
not  instantaneously  set  to  zero,  as  with  the  previous  model.  Instead,  their  gradual 
reduction  is  imposed,  qualitatively  consistent  with  results  from  shear  band  modeling 
reported  in  Wright  and  Walter  (1987).  If  s*r,  s*lz,  sqq,  and  s*z  are  the  four 
components  of  deviatoric  stress  computed  by  EPIC-2  in  such  an  element  for  a  given 
time  step,  then  each  of  these  is  altered  to  the  values  given  by 

Sjj  =  s'j  f(£p)  (17) 

where  for  A efaii  =  0  , 


f(ep)  = 


1  ;  £^€[0,  zponset] 
i0  ;  epe(epnset,  oo) 


(18a) 
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while  for  AeJfljl  >  0  , 


1 

f(ep)  =  1 

0 

The  function  f(ep )  is  drawn  in  Figure  37.  Here  e%nse,  is  the  value  of  equivalent 
plastic  strain  ep  that  existed  in  the  element  at  the  time  step  in  which  the  Zener- 
Hollomon  criterion  was  first  met.  zponset  is  therefore  computed  by  the  code.  On  the 
other  hand,  A efaii  is  a  user-prescribed  material  property.  (e£ni,  +  Aefau)  is  the 
level  of  ep  beyond  which  the  element  no  longer  supports  deviatoric  stresses. 

The  post-onset  conditions  imposed  on  an  element  that  has  failed  according  to 
the  tension  model  are  more  stringent  than  those  imposed  on  elements  that  have 
failed  by  the  shear  band  model.  In  the  former  case  hydrostatic  tensile  stress  as  well 
as  deviatoric  stresses  are  affected,  and  stress  reduction  to  zero  is  instantaneous.  For 
this  reason,  elements  that  have  failed  by  shear  banding  continue  to  be  checked  for 
tensile  failure. 

8.2  Problem  Series  III  Through  VIII.  Six  problem  series  are  run  with  the  tensile 
and  shear  band  failure  models  both  operative.  In  each  series  parameter  pfail  is  fixed 
at  -3.00  GPa  for  both  RHA  and  copper,  and  A efaii  is  fixed  at  either  0.00,  0.10, 

0.25,  0.50,  0.75  or  1.00.  tperode  is  again  varied  over  the  range  of  0.25,  0.50,  0.75, 
1.25,  and  2.00  within  each  series,  with  the  same  value  for  £perode  assigned  to  both 
slidelines. 

Problem  Series  III,  IV,  V,  VI,  VII  and  VIII  correspond  to  a  A efaii  of  0.00, 

0.10,  0.25,  0.50,  0.75  and  1.00,  respectively.  Problem  Series  III,  for  which 
A zpfaii=  0.00,  is  the  case  of  instantaneous  reduction  of  deviatoric  stresses  to  zero 
upon  shear  band  failure  in  an  element.  This  is  the  case  governed  by  Equation  (18a). 
At  the  other  extreme.  Problem  Series  VIII,  for  which  Aefait=  1.00,  allows  gradual 
stress  reduction  following  satisfaction  of  the  Zener-Hollomon  condition,  as 
governed  by  Equation  (18b).  In  this  series  zero  stresses  are  not  imposed  within  an 
element  until  ep  exceeds  zponSet  by  1.00.  Results  for  rhole,  AM  and  t  95  from  the  six 


pp  _  Pp 
S _ ^  onset 

A  £fail 


;  ePefO,  ZPonset\ 

»  onset’  2 onset  Aefaii) 

>  £P  €  [£  onset  +  A  E^an ,  °°) 


(18b) 
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problem  series  are  presented  in  Table  5. 

Figure  38  plots  rhoie  vs.  z%rode  for  the  six  series.  Here  Aefaii  is  a  constant 
parameter  for  each  curve.  On  the  other  hand  Figure  39  plots  rhoie  vs.  Aefaii  with 
£  erode  a  constant  parameter  for  each  curve.  For  each  series,  corresponding  to  a  fixed 
Atpfaii,  rhoie  is  seen  to  increase  with  increasing  £%rode  for  zperode  in  the  range  of  0.25  to 
2.00.  Also,  r foie  decreases  with  increasing  A zfail  for  zperode  in  the  range  of  0.25  to 
2.00,  with  a  single  exception:  at  an  zperode  of  2.00,  rho[e  increases  when  AefaU  is 
increased  from  0.75  to  1.00.  With  this  one  exception  ignored,  the  finding  means 
that  the  more  gradually  deviatoric  stresses  are  reduced  following  the  onset  of  shear 
band  failure,  the  larger  the  final  hole  size. 

Activation  of  the  shear  band  failure  model  allows  for  a  large  increase  in 
computed  hole  size  over  those  obtained  in  Problem  Series  I  and  II.  For  the  series 
with  A efau  equal  to  0.50,  0.75,  and  1.0,  rhole  is  still  consistently  at  least  slightly 
smaller  than  the  experimental  value  of  18.2  mm.  However,  the  curves  in  Figure  38 
corresponding  to  Problem  Series  III,  IV  and  V  and  respective  Azfaii  values  of  0.00, 
0.10,  and  0.25  intersect  the  experimental  result  at  an  zperode  of  about  0.5,  0.7,  and 
1.0,  respectively.  Four  problems  that  produce  reasonably  close  agreement  with 
experimental  rhote ,  arranged  in  order  of  increasing  Azfau,  are  Problem  III-B 
erode =  0.50,  A zfaii=  0.00),  Problem  IV-C  ( zperode=  0.15,  Azfail=  0.10),  Problem 
V-D  ( zper0de=  1.25,  A zfaii-  0.25),  and  Problem  VI-E  (zperode=  2.00,  Azfaii=  0.50). 
These  four  problems  all  produce  rho[e  values  within  4%  of  the  experimental  result. 

Mesh  plots  are  now  used  to  observe  the  evolution  in  time  of  the  solution  to 
Problem  III-B.  Plots  at  1.03,  2.01,  3.01,  4.03,  5.00,  10.01,  20.04,  50.04,  100.06, 
150.05,  200.06,  250.04,  300.06,  450.05  and  500.02  ps  after  impact  are  presented  in 
Figures  40  through  54,  respectively.  The  blue  elements  in  Figure  40  indicate  those 
copper  elements  that  have  failed  in  tension.  The  green  elements  of  RHA  in  Figures 
40  through  54  have  failed  by  shear  banding,  and  the  aquamarine  elements  of  RHA 
in  Figures  42,  43  and  44  have  failed  successively  by  shear  banding  followed  by 
tension.  Figures  40,  41  and  42  show  the  penetration  process  prior  to  perforation. 

In  these  figures  the  target’s  hole  is  surrounded  by  RHA  material  that  has  failed  by 
shear  banding.  By  3.01  ps,  this  failed  region  has  extended  to  the  exit  surface.  A 
single  aquamarine  element  of  RHA  is  located  on  the  exit  surface  in  Figure  42, 
indicating  that  the  phenomenon  of  spallation  has  occurred  slightly  before 
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perforation.  Perhaps  formation  of  the  relatively  small,  early-time  debris  fragments 
noted  in  Section  3  on  the  basis  of  Figure  12  can  therefore  be  attributed  to  the  failure 
mechanism  of  tensile  voids.  According  to  Figures  42  and  43,  perforation  occurs 
between  3.01  and  4.03  (j. s. 

Once  an  element  of  RHA  material  adjacent  to  the  current  hole  boundary  has 
failed  by  shear  banding,  the  stress  reduction  scheme  of  Equations  (17)  and  (18) 
renders  it  prone  to  large  deformation  and  hence  susceptible  to  slideline  erosion. 
Figures  40  through  54  show  that  in  Problem  III-B,  every  RHA  element  that  fails  by 
shear  banding  and  that  does  not  have  at  least  two  of  its  edges  bounded  by  elements 
that  have  not  failed  is  discarded  via  slideline  erosion  by  500  ps.  A  comparison  of 
Figures  48,  49  and  50  reveals  that  a  substantial  amount  of  RHA  material  is  still 
eroding  and  contributing  to  the  final  hole  size  between  100  and  200  ps  after  impact. 
It  is  suggested  that  the  relatively  large,  late-time  fragments  noted  in  Section  3  to  be 
present  in  the  radiographs  of  Figure  12  are  primarily  attributeable  to  the  failure 
mechanism  of  shear  banding  rather  than  tensile  voids.  Figures  51  and  52  show  that 
the  last  failure  by  shear  banding  occurs  between  250  and  300  ps  in  a  single  element 
located  near  the  entrance  surface.  This  element  and  one  of  its  neighbors  erode  from 
the  target  plate  at  between  450  and  500  ps  after  impact,  as  is  clear  from  Figures  53 
and  54.  This  is  the  final  occurrence  of  slideline  erosion  observed  in  Figures  40 
through  54. 

The  mesh  at  500  ps  shown  in  Figure  54  is  taken  to  be  the  final  state  for 
Problem  III-B.  This  figure  therefore  contains  a  prediction  of  the  final  hole 
geometry  and  the  final  spatial  distribution  of  shear  bands  within  a  cross-section. 
Comparison  with  Figure  11  shows  that  this  predicted  hole  geometry  agrees  quite 
closely  with  the  experimental  result,  both  in  terms  of  the  average  hole  radius  rhoie 
and  the  general  shape  of  the  hole  profile.  In  particular  Figure  54  contains  a 
prominent  protuberance  that  juts  into  the  hole  to  form  a  distinguishable  throat  at 
about  the  midsurface  elevation. 

Mesh  plots  at  500  ps  from  Problems  IV-C,  V-D  and  VI-E,  which  also  agree 
closely  with  experiment  in  terms  of  rhoie,  are  shown  in  Figures  55,  56  and  57, 
respectively.  As  with  Problem  III-B,  green  elements  indicate  shear  banding  failure, 
and  the  single  aquamarine  element  in  Figure  57  indicates  shear  banding  failure 
followed  by  tensile  void  failure.  Figures  54  through  57  reveal  a  trend  in  the  number 
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of  RHA  elements  that  have  failed  by  shear  banding  and  that  remain  attached  to  the 
plate  at  500  ps.  As  Aefaii  is  increased  in  the  range  of  0.00  to  0.50,  the  number  of 
these  elements  also  increases.  In  other  words  the  number  of  RHA  elements  able  to 
escape  slideline  erosion  at  500  ps  increases  as  the  rate  at  which  deviatoric  stresses 
are  reduced  following  the  onset  of  shear  banding  decreases.  Note  that  e^rode  is  not 
held  fixed  in  these  four  problems.  zpeTOde  equals  0.50,  0.75,  1.25  and  2.00  in 
Problems  III-B,  IV-C,  V-D  and  VI-E,  respectively.  What  is  held  approximately 
fixed  is  rhoie. 

Next  the  mesh  plots  at  500  ps  from  Problems  III-B,  IV-C,  V-D  and  VI-E  are 
examined  in  terms  of  their  predicted  locations  of  tensile  voids  and  shear  bands. 
These  observations  are  compared  with  the  experimental  results  from  Rd.  10771  in 
Figure  27.  Figures  54,  55  and  56  contain  no  elements  that  have  failed  by  tensile 
voids,  a  result  that  is  clearly  inconsistent  with  the  experimental  presence  of  voids  in 
Figure  27.  Figure  57  contains  only  a  single  element  that  has  failed  in  this  manner, 
so  that  Problem  VI-E  has  also  yielded  an  insufficiently  large  region  of  plate  cross- 
section  that  has  failed  by  voids.  However,  the  location  of  this  single  element,  near 
the  intersection  of  the  hole  boundary  with  the  plate’s  midsurface,  does  constitute  a 
certain  degree  of  agreement  with  experiment,  since  in  Figure  27  all  voids  reside 
within  the  midsurface  protuberance. 

The  calculations  are  more  successful  in  terms  of  predictions  of  shear  band 
locations  in  the  final  cross-section.  Figure  54  from  Problem  III-B  does  display  some 
shear  banding  within  the  midsurface-region  protuberance  that  forms  the  throat. 

This  is  the  region  in  which  most  shear  bands  are  observed  in  Figure  27.  However, 
Figure  54  does  not  display  shear  banding  along  the  boundary  of  this  protuberance, 
where  it  is  conspicuous  in  Figure  27.  Figure  54  also  displays  shear  banding  near  the 
intersection  of  the  hole  boundary  and  the  entrance  surface,  consistent  with  Figure 
27.  Two  elements  of  shear  banding  are  also  contained  in  Figure  54  near  the 
intersection  of  the  hole  boundary  and  the  exit  surface.  This  last  result  is 
inconsistent  with  Figure  27.  However,  the  caveat  stated  in  Section  6  should  be 
recalled.  The  shear  bands  denoted  in  Figure  27  are  those  that  were  observable  with 
an  optical  microscope  at  1000X  magnification.  Perhaps  additional  shear  bands  could 
have  been  observed  in  these  and  other  regions  of  the  cross-section  by  more 
discriminating  means.  Thus  in  effect  Figure  27  constitutes  a  lower  bound  on  the 
number  of  shear  bands  present  in  the  experiment. 
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Figure  55  from  Problem  IV-C  contains  a  cylindrical  ring  about  5  mm  thick  in 
which  shear  bands  are  distributed  rather  uniformly.  Thus  there  is  little  dependence 
on  the  density  of  shear  bands  with  elevation  z  predicted  for  this  problem.  This  is 
inconsistent  with  the  preponderance  of  shear  bands  located  near  the  midsurface 
elevation  in  Figure  27.  The  number  of  shear-banded  elements  within  the  cylindrical 
ring  is  about  four  times  the  total  number  contained  in  Figure  54.  The  general 
impression  is  that  a  larger  portion  of  the  cross-section  in  Figure  55  has  failed  by 
shear  banding  than  in  the  experimental  cross-section  of  Figure  27,  but  the  caveat  of 
Section  6  must  again  be  borne  in  mind.  Figure  55  contains  a  single  shear-banded 
element  at  the  exit  surface  and  about  18  mm  from  the  hole  boundary.  Figure  56 
from  Problem  V-D  shows  a  hole  profile  that  is  "flatter",  or  has  less  z  dependence, 
than  in  Figures  54  and  55.  Shear  banding  is  predicted  in  Figure  56  along  most  of 
the  hole  boundary.  Within  a  5-mm-thick  cylindrical  ring  adjacent  to  the  hole 
boundary,  a  larger  percentage  of  elements  have  shear  banded  in  Figure  56  than  in 
Figure  55.  In  Figure  57  from  Problem  VI-E,  the  corresponding  5-mm-thick  ring  has 
almost  entirely  failed  by  shear  banding.  In  both  Figures  56  and  57,  the  portion  of 
the  cross-section  that  has  failed  by  shear  banding  clearly  exceeds  the  experimental 
result  of  Figure  27. 

Figures  58  through  63  plot  cumulative  mass  lost  by  the  target  plate  as  a  function 
of  time  after  impact  for  Problem  Series  III,  IV,  V,  VI,  VII  and  VIII,  respectively. 

In  Figure  58  the  curve  for  which  zperode  =  0-50  corresponds  to  Problem  III-B.  This 
curve  can  be  divided  into  four  regions.  For  the  first  25  ps  after  impact,  mass  loss 
occurs  at  the  rapid  average  rate  of  1.9  g  ps'1.  The  average  mass  loss  rate  between 
25  and  100  ps  is  0.52  g  ps1.  From  100  to  200  ps,  the  average  rate  further  reduces 
to  the  still  appreciable  0.17  g  ps1.  This  is  consistent  with  the  continuing  mass  loss 
detected  in  the  mesh  plots  of  Figures  48,  49  and  50.  By  200  ps  the  curve  has 
reached  a  plateau;  no  additional  erosion  occurs  until  489.8  ps,  consistent  with  the 
observations  based  on  the  mesh  plots  in  Figures  50  through  54. 

The  value  attained  at  500  ps  by  each  curve  in  Figures  58  through  63  is 
identified  with  AM,  the  cumulative  mass  that  the  plate  has  lost  during  the  500  ps 
following  initial  impact.  These  results  for  AM  are  plotted  vs.  zperode  with  Aefau  a 
parameter  in  Figure  64  and  vs.  Ae^0 17  with  zperodr  a  parameter  in  Figure  65. 
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Figures  64  and  65  show  that  in  each  of  the  Problem  Series  III  through  VIII,  AM 
increases  with  increasing  zperode  for  zperode  in  the  range  of  0.25  to  2.00.  These  figures 
also  show  that  for  eprode  fixed  at  0.25,  0.50,  0.75,  1.25  or  2.00,  AM  decreases  as 
Ae%u  is  increased  in  the  range  of  0.00  to  1.00.  AM  has  been  dramatically  increased 
over  the  values  obtained  in  Problem  Series  I  and  II  by  the  activation  of  the  shear 
band  model.  Thus  whereas  AM  results  from  Problem  Series  I  and  II  never  exceeded 
28.6  g,  in  Figure  64  the  curves  pertaining  to  Problem  Series  III,  IV  and  V  intersect 
the  experimental  result  of  136  g  at  an  zperode  of  about  0.8,  1.1  and  1.6,  respectively. 
The  curves  for  Problem  Series  VI,  VII  and  VIII  remain  below  this  experimental 
result  throughout  the  strode  range  of  0.00  to  2.00. 

Recall  that  similar  observations  were  made  for  rho[e  on  the  basis  of  Figures  38 
and  39,  which  resemble  Figures  64  and  65.  In  Figure  38  the  curves  of  rhoie  vs.  zperode 
from  Problem  Series  III,  IV  and  V  also  intersected  the  experimental  result  of  18.2 
mm,  but  at  the  lower  zperode  values  of  about  0.5,  0.7  and  1.0,  respectively.  As  a 
consequence  of  this  difference,  Problems  III-B,  IV-C,  V-D  and  VI-E,  which  all 
produced  rhote  values  within  4%  of  the  experimental  result,  are  now  seen  in  Table  5 
to  produce  AM  values  that  are  smaller  than  the  experimental  result  by  22%,  18%, 
18%  and  30%,  respectively.  These  discrepancies  between  computed  AM  values 
from  Problems  III-B,  IV-C,  V-D  and  VI-E  and  the  experimental  result  of  136  g  are 
attributed  to  the  removal  of  insufficient  RHA  material  in  these  four  problems, 
despite  the  fact  that  the  final  hole  radius  displays  close  agreement  with  experiment 
in  each  case.  Apparently  too  much  RHA  material  has  been  compressed  or  bent  and 
retained  in  the  calculations  instead  of  being  discarded  by  slideline  erosion. 

The  problem  which  produces  a  AM  result  in  closest  agreement  with  experiment 
is  Problem  III-C.  AM  for  this  problem  is  130.3  g  (Table  5),  which  is  only  4%  less 
than  the  experimental  result  of  136  g.  The  rhole  value  obtained  from  Problem 
III-C  is  20.2  mm,  which  is  11%  larger  than  the  experimental  result  of  18.2  mm. 
Thus  the  correct  amount  of  RHA  material  has  been  retained  in  Problem  III-C,  but 
the  final  hole  radius  is  somewhat  larger  than  in  the  experiment. 

The  four  problems  that  produce  AM  results  in  closest  agreement  with 
experiment  after  Problem  III-C  are  Problems  IV-C,  IV-D,  V-D  and  V-E.  The  AM 
values  computed  for  these  problems  agree  with  the  experimental  measurement  from 
Rd.  4190  to  within  18%,  12%,  18%  and  16%,  respectively.  However,  as  with 
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Problem  III-C,  these  four  problems  also  produce  rhole  values  larger  than  the 
experimental  result  of  Equation  (1),  by  3%,  20%,  4%  and  23%,  respectively. 

Recall  that  Problems  IV-C  and  V-D  were  previously  noted  to  have  produced  rhoU 
values  in  good  agreement  with  experiment.  The  meshes  at  500  ps  for  these  two 
problems  were  presented  in  Figures  55  and  56. 

Return  now  to  the  mass  loss  vs.  time  plots  of  Figures  58  through  63,  applicable 
to  Problem  Series  III  through  VIII,  respectively.  When  these  are  compared  with  the 
mass  loss  vs.  time  plots  from  Problem  Series  I  and  II,  contained  in  Figures  23  and 
36,  respectively,  it  becomes  apparent  that  activation  of  the  shear  band  model  has 
greatly  extended  the  time  at  which  target  erosion  comes  to  completion.  In  fact  on 
the  basis  of  Figures  58  through  63,  some  of  the  problems  in  Problem  Series  III 
through  VIII  still  appear  to  be  undergoing  gradual  additional  erosion  of  RHA  at  500 
ps  after  impact.  Results  for  r  95,  the  time  at  which  the  plate  has  lost  95%  of  mass 
AM,  are  listed  in  Table  5.  All  t .95  listings  for  Problem  Series  III  through  VIII 
satisfy  the  experimental  bound  of  Equation  (3).  t  95  results  in  Table  5  do  not 
exhibit  a  simple  dependence  on  either  zper0de  or 

9.  EXAMINATION  OF  RESULTS’  SENSITIVITY  TO  MESH  FINENESS 

All  calculations  thus  far  have  employed  the  initial  meshes  for  the  jet  particle 
and  the  RHA  target  plate  that  are  shown  in  Figure  15.  A  question  that  arises  is  to 
what  extent  the  results  presented  above  for  the  eight  problem  series  have  been 
spuriously  determined  by  the  geometry  of  this  mesh,  rather  than  by  the  modeling  of 
the  various  physical  phenomena  in  the  problem.  In  order  to  pursue  this  question. 
Problems  I-C,  II-C,  and  III-B  will  be  rerun,  this  time  using  the  same  mesh  for  the 
jet  particle  as  in  Figure  15,  but  with  a  finer  RHA  target  mesh.  The  new  mesh, 
shown  in  Figure  66,  consists  of  isosceles  right  triangles  that  are  similar  to  those  of 
Figure  15,  but  with  their  longest  edge  length  equal  to  0.65  mm,  or  one-half  the 
longest  edge  length  of  each  target  element  in  Figure  15.  Thus  the  total  number  of 
RHA  elements  has  been  quadrupled  over  that  in  Figure  15,  from  3200  to  12800. 
There  are  now  twenty  layers  of  crossed  triangles  through  the  plate’s  thickness. 

These  new  problems,  employing  the  mesh  of  Figure  66,  are  called  Problems  I-C(R), 
II-C(R),  and  III-B(R). 
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Problem  I-C(R)  has  neither  failure  model  activated  and  cutoff  erosion  strain 
£erode  set  to  0.75  for  both  slidelines.  Mesh  plots  from  this  problem  at  3.00,  4.00 
and  500.02  ps  after  impact  are  displayed  in  Figures  67,  68  and  69,  respectively. 
These  can  be  compared  with  Figures  16,  17  and  18,  which  apply  to  Problem  I-C  at 
similar  times.  Hole  geometries  from  the  two  problems  are  in  reasonable  agreement 
at  the  two  earlier  times.  Perforation  occurs  at  between  3  and  4  ps  in  both  cases.  At 
500  ps,  the  final  average  hole  radius,  rhote,  is  slightly  smaller  in  Problem  I-C(R) 
than  in  Problem  I-C,  and  hence  a  little  further  from  the  experimental  result  of  18.2 
mm.  rhoie  is  determined  to  be  9.7  mm  from  Figure  69  and  11.0  mm  from  Figure  18, 
for  a  difference  of  13%.  In  addition  the  mesh  at  500  ps  contains  more  prominent 
lips  at  the  intersections  of  the  hole  boundary  with  the  entrance  and  exit  surfaces  in 
Problem  I-C(R)  than  in  Problem  I-C. 

The  total  mass  lost  by  the  target  plate  is  plotted  as  a  function  of  time  for  both 
Problems  I-C  and  I-C(R)  in  Figure  70.  After  the  initial  3  ps,  this  quantity  is  seen  to 
be  consistently  smaller  in  Problem  I-C(R)  than  in  Problem  I-C.  The  total  mass  loss 
at  500  ps,  AM,  is  21.3  g  in  the  former  and  27.7  g  in  the  latter  problem.  The 
difference  here  is  30%.  The  time  at  which  95%  of  AM  has  been  lost,  t  95,  is 
determined  to  be  49.9  ps  in  Problem  I-C(R).  This  is  210. %  larger  than  the  value  of 
16.1  ps  determined  for  Problem  I-C,  but  still  smaller  than  the  experimental  bound 
of  Equation  (3).  Results  for  rhole,  AM  and  t  95  from  Problem  I-C(R)  are  collected 
in  Table  6. 

Problem  II-C(R)  has  the  tensile  void  failure  model  activated  for  OFHC  copper 
and  RHA  steel,  with  material  parameter  pfaii  set  to  -3.00  GPa  for  both  metals.  The 
cutoff  strain  for  slideline  erosion,  zper0de >  is  set  to  0.75  for  both  slidelnes,  as  in 
Problem  I-C(R).  Mesh  plots  from  Problem  II-C(R)  at  1.00,  3.01  and  500.00  ps 
after  impact  are  presented  in  Figures  71,  72  and  73,  respectively.  These  are 
comparable  with  Figures  28,  29  and  30,  showing  mesh  plots  from  Problem  II-C  at 
corresponding  times.  In  general  good  qualitative  agreement  with  the  earlier 
problem  is  obtained.  At  1.00  ps,  failed  copper  elements  are  seen  near  the  leading 
edge  of  the  jet  particle.  In  contrast  to  Problem  II-C,  a  single  failed  element  of 
target  material  occurs  as  well,  at  the  entrance  surface  of  the  hole  boundary.  At  3.01 
I  s,  the  phenomenon  of  spallation,  as  visualized  by  the  presence  of  failed  target 
elements  near  the  exit  surface,  is  observed  again.  The  spall  fracture  surface  seems 
to  be  recessed  by  very  roughly  one  element  layer,  or  0.6  mm,  from  the  exit  surface 
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in  Problem  II-C(R),  an  effect  that  could  not  be  resolved  with  the  coarser  mesh  of 
Problem  II-C.  At  500  (is,  both  Figures  30  and  73  reveal  a  small  region  of  RHA 
that  has  failed  in  tension.  This  region  is  located  near  the  intersection  of  the  hole 
boundary  with  the  plate’s  midsurface  in  both  figures.  As  mentioned  in  Section  7, 
this  location  of  tensile  failure  is  consistent  in  a  general  sense  with  Figure  27,  which 
shows  voids  scattered  about  the  midsurface  protuberance.  The  hole  radius  at  500  ps 
in  Problem  II-C(R)  agrees  closely  with  that  from  Problem  II-C  from  the  exit  surface 
to  the  vicinity  of  the  midsurface.  Between  the  midsurface  and  the  entrance  surface, 
the  hole  radius  at  500  ps  is  slighly  smaller  in  Problem  II-C(R)  than  in  Problem  II-C. 
rhoie,  the  radius  averaged  over  the  thickness,  is  10.2  mm  in  Problem  II-C(R),  as 
opposed  to  10.9  mm  in  Problem  II-C.  The  difference  between  rhole  values  from 
Problems  II-C(R)  and  II-C  is  7%.  Both  are  much  smaller  than  the  experimental 
value  of  18.2  mm. 

A  comparison  of  Figures  69  and  73  shows  that  with  this  finer  mesh,  activation 
of  the  tensile  void  failure  model  has  caused  the  final  hole  radius  to  grow  slightly 
throughout  most  of  the  target  plate  thickness.  The  average  radius,  rhole,  is  9.7  mm 
in  Problem  I-C(R)  and  10.2  mm  in  Problem  II-C(R).  Recall  from  Section  7  that 
with  the  coarser  mesh,  a  very  small  decrease  in  rhoie  was  found  to  accompany 
activation  of  this  model;  rho[e  was  11.0  mm  in  Problem  I-C  and  10.9  mm  in  Problem 
II-C. 


Figure  74  presents  plots  of  cumulative  mass  lost  by  the  target  plate  vs.  time  for 
Problems  II-C  and  II-C(R).  According  to  this  figure,  the  mass  loss  is  smaller  in 
Problem  II-C(R)  than  in  Problem  II-C  at  all  times  after  5  ps.  This  is  reminiscent 
qualitatively  of  the  relationship  between  mass  loss  results  from  Problems  I-C  and 
I-C(R).  Quantitatively,  however,  the  difference  between  mass  loss  results  obtained 
with  the  finer  mesh  and  with  the  coarser  mesh  is  less  in  Figure  74  than  in  Figure  70 
throughout  this  time  range.  The  cumulative  mass  lost  at  500  ps,  AM,  is  seen  in 
Figure  74  to  be  27.2  g  for  Problem  II-C  and  23.1  g  for  Problem  II-C(R),  for  a 
difference  of  18%.  Both  are  much  smaller  than  the  experimental  result  of  136  g. 
t  95,  the  time  at  which  the  target  plate  has  lost  cumulative  mass  equal  to  95%  of  AM, 
equals  22.7  and  28.1  ps  for  Problems  II-C  and  II-C(R),  respectively.  These  times 
differ  from  each  other  by  24%  and  are  both  much  earlier  than  the  experimental 
bound  of  Equation  (3).  Note  that  AM  from  Problem  II-C(R)  is  slightly  larger  than 
the  value  obtained  from  Problem  I-C(R).  This  is  consistent  with  the  relationship 
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observed  for  rho[e  from  the  two  problems,  t  95  from  Problem  II-C(R)  is  less  than 
the  result  from  Problem  I-C(R),  and  both  are  substantially  smaller  than  the 
experimental  bound.  Results  for  r^0/e,  AM  and  r.95  obtained  from  Problem  II-C(R) 
are  added  to  Table  6. 

Problem  III-B(R)  has  the  shear  band  failure  model  activated  for  the  RHA  in 
addition  to  the  tensile  void  failure  model  for  both  OFHC  copper  and  RHA  with  pfait 
set  to  -3.00  GPa.  Shear  band  parameter  Aefaii  is  set  to  0.00,  and  cutoff  erosion 
strain  e^r0tJe  is  set  to  0.50  for  both  slidelines.  Recall  that  Problem  III-B,  the 
counterpart  of  Problem  III-B(R)  with  the  coarser  mesh,  gave  particularly  satisfying 
agreement  with  experiment  in  terms  of  the  final  hole  geometry.  A  question  that 
comes  up  now  is  to  what  extent  that  agreement  is  preserved  with  the  finer  mesh  in 
Problem  III-B(R).  Mesh  plots  from  Problem  III-B(R)  are  presented  at  1.01,  2.01, 
3.01,  4.00,  5.02,  10.02,  20.01,  50.03,  100.01,  200.02  and  500.01  ps  after  impact  in 
Figures  75  through  85.  These  can  be  compared  with  Figures  40  through  48,  50  and 
54,  which  pertain  to  similar  times  for  Problem  III-B.  Mesh  plots  at  times  prior  to 
perforation  are  in  good  qualitative  agreement.  Spallation  in  the  target  plate  is  again 
observed  at  3  ps  for  both  meshes.  The  radial  extent  of  the  spall  fracture  surface  is 
about  4  mm  in  Problem  III-B(R)  and  only  about  1  mm  in  Problem  III-B.  Also,  the 
spall  fracture  surface  is  recessed  about  0.8  mm  from  the  exit  surface  in  Problem 
III-B(R),  while  no  such  recession  can  be  resolved  with  the  coarser  mesh  of  Problem 
III-B. 

Following  perforation,  mesh  plots  from  Problem  III-B(R)  continue  to  display 
good  qualitative  agreement  with  those  from  Problem  III-B  through  10  ps,  both  in 
terms  of  the  hole  boundary’s  geometry  and  the  general  location  of  zones  of  elements 
that  have  failed.  At  4  and  5  ps,  a  small  zone  of  elements  that  have  failed  in  tension 
subsequent  to  failure  by  shear  banding  exists  at  roughly  the  same  location  in  both 
problems.  Furthermore,  by  10  ps  this  zone  has  eroded  away  in  both  problems. 

By  20  ps  after  perforation,  differences  in  hole  profiles  from  Problems  III-B 
and  III-B(R)  become  apparent  and  seem  to  be  driven  by  differences  in  the  spatial 
distributions  of  zones  of  RHA  that  have  failed  by  shear  banding.  For  example,  as 
was  noted  in  Section  8,  a  midsurface-level  protuberance  remains  present  in  the  mesh 
at  500  ps  in  Problem  III-B  (Figure  54).  This  feature  is  in  good  agreement  with  the 
experimental  hole  profile  of  Figure  11  and  shows  promise  of  being  preserved  in 
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Problem  III-B(R)  on  the  basis  of  the  mesh  plot  at  20  (is  (Figure  81).  In  this  figure 
a  group  of  RHA  elements  which  have  not  failed  is  seen  to  jut  out  into  a  region  of 
failed  elements  which  line  the  hole  boundary.  However,  erosion  of  this  group  of 
unfailed  elements  occurs  between  50  and  100  ps  (Figures  82  and  83).  This  erosion 
is  assisted  by  the  presence  of  contiguous  failed  elements  located  at  about  midsurface 
elevation  and  radially  outward  from  the  group  of  unfailed  elements  in  question. 

The  result  is  that  no  prominent  midsurface  protuberance  is  present  in  the  final  mesh 
from  Problem  III-B(R)  (Figure  85).  Instead,  the  hole’s  main  throat,  or  most 
narrow  cross-section,  occurs  near  the  entrance  surface,  with  a  lesser  throat  occurring 
near  the  exit  surface.  The  hole  is  relatively  wide  at  midsurface  elevation.  rhoie,  the 
hole’s  radius  averaged  over  the  plate  thickness  at  500  ps,  is  determined  for  Problem 
III-B(R)  on  the  basis  of  nodal  coordinates  along  the  red  portion  of  the  target 
boundary  in  Figure  86.  The  result  is  17.3  mm  for  rho[e.  This  agrees  reasonably  well 
with  the  experimental  result  of  18.2  mm,  despite  the  large  difference  between  hole 
profile  shapes  in  Figure  85  or  86  and  Figure  11.  It  also  agrees  reasonably  well  with 
the  result  of  18.1  mm  from  Problem  III-B.  The  difference  amounts  to  a  5% 
reduction  in  rho[e  corresponding  to  a  doubling  in  mesh  fineness. 

Next  the  locations  of  failed  elements  in  Figure  85  from  Problem  III-B(R)  will 
be  examined  in  light  of  the  mesh  in  Figure  54  from  Problem  III-B  and  of 
experimental  results  from  Rd.  10771  that  were  discussed  in  Section  6.  First  note 
from  Figure  85  that  no  elements  that  have  failed  by  tensile  voids  have  escaped 
slideline  erosion  by  500  ps.  This  is  inconsistent  with  the  experimental  result  of 
Figure  27,  which  shows  at  least  nine  large,  discrete  voids  distributed  around  the 
midsurface  protuberance.  Moreover  there  are  no  elements  that  have  failed  by 
tensile  voids  in  the  mesh  at  50  ps,  either  among  elements  in  the  midsurface  that  will 
shortly  erode  away  or  elsewhere  (Figure  82).  On  the  other  hand  the  mesh  in  Figure 
85  has  retained  numerous  RHA  elements  that  have  failed  by  shear  banding.  These 
occur  mainly  in  the  vicinity  of  the  hole  boundary  and  between  the  entrance  surface 
and  the  midsurface.  In  contrast  Figure  54  shows  a  more  uniform  distribution  of 
predicted  shear  bands  with  respect  to  elevation  in  the  plate.  Figure  27  from  the 
experiment  shows  two  or  three  shear  bands  detected  near  the  entrance  surface, 
consistent  with  a  cluster  of  failed  elements  in  Figure  85.  Most  of  the  shear  bands  in 
Figure  27  lie  within  the  midsurface  protuberance.  This  feature  of  the  experimental 
distribution  is  more  consistent  with  the  mesh  from  Problem  III-B(R)  at  50  ps  than 
at  500  p  s. 
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Mass  lost  by  the  target  plate  is  plotted  as  a  function  of  time  in  Figure  87  for 
Problems  III-B  and  III-B(R).  At  all  times  after  5  ps,  the  mass  loss  in  Problem 
III-B(R)  is  less  than  that  in  Problem  III-B.  At  500  ps,  the  net  mass  loss,  AM,  is 
97.2  g  in  Problem  III-B(R)  and  105.5  g  in  Problem  III-B  (Table  5),  corresponding 
to  an  8.5%  difference.  This  difference  of  8.3  g  bounds  the  difference  at  all  earlier 
times.  The  percentage  difference  between  AM  from  the  two  problems  is  larger  than 
that  found  above  for  the  quantity  rhoie.  This  last  comparison  is  correlated  with  the 
flatter  hole  profile  (i.e.,  less  z  dependence  of  hole  radius)  at  500  ps  for  Problem 
III-B(R)  than  for  Problem  III-B. 

t. 95  is  found  from  Figure  87  to  be  158.3  ps  for  Problem  III-B(R).  This 
represents  a  13.7%  decrease  from  the  result  of  183.4  ps  from  Problem  III-B  (Table 
5).  These  r  95  results  from  both  problems  satisfy  the  experimental  bound  of 
Equation  (3). 

The  results  for  rhoU,  AM  and  t  95  from  Problem  III-B(R)  are  added  to  Table  6. 
Activation  of  the  shear  band  failure  model  with  the  finer  mesh  has  produced  some 
of  the  same  general  effects  as  in  Section  8  with  the  coarser  mesh:  The  final 
perforation  hole  has  grown  substantially  larger  and  has  been  placed  in  closer 
agreement  with  experiment  in  terms  of  size,  as  witnessed  by  results  for  rhote  and 
AM.  Also,  the  duration  of  the  hole  growth  process,  as  measured  by  t .95,  has  been 
significantly  increased  over  the  results  obtained  with  neither  failure  model  or  with 
only  the  tensile  void  model  active  and  has  been  placed  in  agreement  with  the  bound 
derived  from  the  radiographs  of  Figure  12.  The  shape  of  the  final  hole  profile  from 
Problem  III-B(R)  is  flatter  than  both  the  results  from  Problem  III-B  and  from  Rd. 
10771.  However,  variation  of  the  parameters  zperode  and  Atfaii  in  series  of  problems 
with  the  finer  mesh  could  presumably  have  recovered  agreement  with  the 
experimental  hole  shape. 

The  final  issue  to  consider  is  the  added  CPU  times  associated  with  the  finer 
mesh.  All  calculations  described  in  this  report  have  been  performed  on  the  Cray 
X-MP/48  computer  operated  by  the  U.S.  Army  Ballistic  Research  Laboratory.  The 
CPU  times  expended  in  running  Problems  I-C,  I-C(R),  II-C,  II-C(R),  III-B  and 
III-B(R)  out  to  500  ps  are  listed  in  Table  7.  CPU  times  for  Problems  I-C  and  II-C 
are  very  similar,  as  are  those  for  Problems  I-C (R )  and  II-C(R).  Thus  for  both  the 
coarser  and  the  finer  target  mesh,  activation  of  only  the  tensile  void  failure  model 
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has  been  accompanied  not  only  by  little  effect  on  results  for  rhoie,  AM  and  t  95,  but 
also  by  little  increase  in  required  CPU  time.  In  contrast,  for  both  meshes  activation 
of  the  shear  band  failure  model  has  been  accompanied  by  substantial  reduction  in 
required  CPU  time.  The  main  reason  for  this  reduction  is  probably  related  to  the 
time  step,  determined  in  EPIC-2  by  means  of  the  Courant  stabililty  condition 
(Johnson  1976).  According  to  this  condition,  at  any  given  time  in  the  calculation  the 
time  step  is  limited  by  the  current  dimensions  of  the  smallest  non-eroded  element. 

In  problems  with  the  shear  band  failure  model  active,  more  target  elements  become 
eroded  and  fewer  highly  distorted  elements  remain  attached  to  the  plate.  This 
prevents  the  time  step  from  becoming  as  small  as  in  the  other  problems.  Doubling 
the  mesh  fineness  is  seen  in  Table  7  to  have  increased  the  required  CPU  time  by 
factors  of  7.7,  7.8  and  9.0  in  the  case  of  Problems  I-C,  II-C  and  III-B,  respectively. 
A  factor  of  8.0  increase  was  anticipated  on  the  basis  of  two  considerations:  (i) 
Halving  the  initial  element  dimensions  very  roughly  halves  the  time  step  throughout 
the  course  of  the  calculation,  (ii)  Quadrupling  the  number  of  elements  very  roughly 
quadruples  the  number  of  operations  required  at  each  time  step. 

The  salient  results  from  this  section  can  be  briefly  summarized  as  follows:  For 
each  of  the  three  problems  run  with  the  finer  mesh  of  Figure  66,  the  final  hole  size, 
as  measured  both  by  rhoU  and  AM,  is  smaller  than  the  result  from  the  counterpart 
problem  with  the  coarser  mesh.  The  differences  in  terms  of  rhole  are  13%,  7%  and 
5%  for  Problems  I-C(R),  II-C(R)  and  III-B(R),  respectively.  The  differences  in 
terms  of  AM  for  these  respective  problems  are  30%,  18%  and  8.5%.  rhoie  and  AM 
results  from  Problems  I-C,  II-C,  I-C(R)  and  II-C(R)  are  all  substantially  smaller 
than  the  experimental  results,  while  results  for  these  quantities  from  Problems  III-B 
and  III-B(R)  are  all  in  reasonable  agreement  with  experiment.  Specifically,  rhou  and 
AM  are  respectively  0.5%  and  30.%  smaller  than  the  experimental  results  in  the  case 
of  Problem  III-B  and  respectively  5%  and  40.%  smaller  in  the  case  of  Problem 
III-B(R).  The  final  hole  shape  from  Problem  III-B  is  also  in  reasonable  agreement 
with  experiment,  but  the  final  hole  profile  from  Problem  III-B(R)  is  flatter  than 
according  to  experiment.  The  quantity  r.95  has  been  increased  in  the  case  of 
Problems  I-C(R)  and  II-C(R)  by  210.%  and  24%,  respectively,  but  decreased  by 
13.7%  in  the  case  of  Problem  III-B(R)  relative  to  the  values  obtained  with  the 
coarser  mesh.  The  results  for  t ,95  from  Problems  I-C,  II-C,  I-C(R)  and  II-C(R)  are 
smaller  than  in  the  experiment,  while  those  from  Problems  III-B  and  III-B(R)  are 
consistent  with  experiment.  Corresponding  to  the  doubling  of  target  mesh  fineness. 
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CPU  times  on  a  Cray  X-MP/48  have  increased  by  factors  of  7.7,  7.8  and  9.0  in  the 
case  of  Problems  I-C,  II-C  and  III-B,  respectively. 

The  above-cited  quantitative  differences  between  results  for  rhote,  A M  and  t. 95 
from  Problems  III-B  and  III-B(R)  provide  measures  of  the  degree  to  which  the 
mesh  of  Figure  15  is  adequate  for  studying  parametric  dependencies  of  the  solution 
upon  £■  erode>  Pfaii  and  A efau,  as  was  done  in  Section  8.  This  degree  of  adequacy  must 
be  weighed  against  the  added  expense  in  CPU  time  associated  with  performing 
parametric  studies  with  the  finer  mesh  of  Figure  66.  For  most  researchers,  CPU 
times  given  in  Table  7  for  problems  involving  the  refined  mesh  are  prohibitive  for 
performing  the  numerous  runs  required  for  parametric  studies.  A  rational  basis  for 
the  use  of  the  mesh  in  Figure  15  is  therefore  indicated. 

10.  CONCLUDING  REMARKS 

This  report  has  discussed  the  modeling  of  a  particular  experiment,  in  which  a 
shaped  charge  warhead  was  fired  at  15.23-charge-diameter  standoff  into  a  13.0-mm- 
thick  RHA  plate.  Two  quantities  evaluated  from  the  experiment  and  that  served  as 
benchmarks  for  the  calculations  are  rhoie,  the  final  hole  radius  averaged  over  the 
plate  thickness,  and  AM,  the  total  mass  lost  by  the  plate.  A  third  quantity,  which 
was  bounded  on  the  basis  of  results  from  two  experiments  quite  similar  to  Rd. 

10771,  is  t  95,  the  time  after  impact  at  which  the  plate  has  lost  95%  of  mass  AM. 

The  experimental  results  obtained  for  rh0[e,  AM  and  t_95  were  18.2  mm,  136  g,  and 
the  bound  t  g5  >  73.9  (xs,  respectively.  A  fourth  benchmark  was  provided  by  a 
diagram  showing  locations  of  prominent  voids  and  shear  bands  in  a  cross-section  of 
the  perforated  plate  from  the  original  experiment. 

A  five-parameter  geometric  model  for  the  leading  jet  particle  was  presented. 

Its  parameters  were  evaluated  on  the  basis  of  a  radiograph.  The  particle’s  length, 
volume  and  centroid  location  were  determined  in  the  process. 

Algorithms  were  developed  to  simulate  material  damage  associated  with  the  two 
micromechanical  phenomena  of  voids  and  shear  bands,  both  observed  in 
micrographs  of  RHA  specimens  from  the  perforated  plate.  Each  algorithm,  or 
failure  model,  consists  of  an  onset  criterion  and  a  post-onset  prescription  on 
stresses.  The  tensile  void  failure  model,  applied  to  both  the  OFHC  copper  of  the 
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leading  jet  particle  and  the  RHA  steel  of  the  target  plate,  introduces  the  parameter 
Pfaii  for  each  of  the  two  metals.  pfau  is  the  negative  value  of  pressure  which  when 
reached  in  an  element  causes  that  element  to  thereafter  be  restricted  from  bearing 
any  deviatoric  stress  or  hydrostatic  tension.  The  shear  band  failure  model,  applied 
only  to  the  steel,  employs  as  its  onset  criterion  the  Zener-Hollomon  condition.  This 
model  introduces  the  single  parameter  Ae^,,/,  the  additional  equivalent  plastic  strain 
which  a  failed  element  can  support  beyond  the  level  it  had  attained  at  the  time  of 
shear  banding  onset.  These  algorithms  were  installed  into  Lagrangian  hydrocode 
EPIC-2. 

A  finite  element  mesh  was  developed  for  the  leading  jet  particle  and  for  the 
target  plate.  That  for  the  target  plate  consisted  of  3200  isosceles  right  triangles, 
each  with  its  longest  edge  equal  to  1.3  mm  in  length.  The  contact  surface  between 
the  meshes  for  the  leading  copper  jet  particle  and  the  RHA  target  plate  were 
modeled  by  means  of  the  eroding  slideline  algorithm  contained  in  EPIC-2.  This 
eroding  slideline  algorithm  introduces  the  material  parameter  e£rode,  a  critical 
equivalent  plastic  strain,  for  each  of  OFHC  copper  and  RHA  steel.  With  these 
initial  meshes  fixed,  eight  series  of  five  problems  each  were  run.  In  each  problem 
Eerode  was  assigned  the  same  value  for  both  materials. 

The  first  series  had  neither  failure  model  active  for  the  two  metals.  The  only 
parameter  varied  was  e£rode,  over  the  range  of  values  0.25,  0.50,  0.75,  1.25  and 
2.00.  All  five  problems  predicted  a  plate  hole  that  is  smaller  and  whose  formation 
is  completed  earlier  than  in  the  experiment.  The  largest  values  obtained  for  rhole, 
AM  and  t  95  in  this  series  were  13.9  mm,  27.7  g,  and  19.4  ps,  respectively, 
corresponding  to  respective  zper0de  values  of  2.00,  0.75,  and  2.00. 

The  experiment  was  then  remodeled  using  EPIC-2  with  only  the  tensile  void 
failure  model  activated.  Throughout  this  second  problem  series,  parameter  pfail  was 
assigned  the  value  -3.00  GPa  for  both  RHA  and  copper.  Activation  of  the  tensile 
void  failure  model  alone  did  not  significantly  increase  rhole  or  AM.  Neither  was  r. 95 
enlarged  significantly,  except  for  the  case  of  e£rorfe  equal  to  2.00. 

Activation  of  both  the  tensile  void  and  the  shear  band  failure  models  for 
Problem  Series  III  through  VIII  did  allow  for  substantial  increases  in  rhoiey  AM  and 
r .95 .  Pfaii  was  again  fixed  at  -3.00  GPa  for  both  metals.  Each  of  the  four  problems 
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for  which  (  zperode,  A tfau)  was  assigned  the  values  (0.50,  0.00),  (0.75,  0.10),  (1.25, 
0.25)  or  (2.00,  0.50)  produced  an  rhoie  value  within  4%  of  the  experimental  result. 
Of  these  four  problems,  the  one  involving  (0.50,  0.00)  produced  the  closest 
agreement  with  experiment  in  terms  of  the  final  hole  profile’s  general  shape  and  the 
final  location  of  shear  bands  within  a  cross-section  of  the  target  plate.  None  of 
these  problems  predicted  a  sufficient  portion  of  the  plate’s  final  cross-section  to  be 
taken  up  by  voids.  These  four  problems  all  yielded  AM  results  that  were  lower  than 
that  of  the  experiment,  by  22%,  18%,  18%  and  30%,  respectively.  Evidently  more 
highly  compressed  and  bent  RHA  material  remained  attached  to  the  plate  in  the 
calculations  than  in  the  experiment.  On  the  other  hand  it  was  found  possible  to 
produce  closer  agreement  with  experiment  in  terms  of  AM  by  tolerating  more 
deviation  in  rhote.  Thus  a  fifth  problem,  for  which  (zperode,  Aefau)  was  set  to  (0.75, 
0.00),  produced  a  AM  result  within  4%  of  experiment,  but  an  rhoie  result  that  was 
11%  larger  than  according  to  experiment.  All  calculations  performed  with  both 
failure  models  active  yielded  t  95  results  consistent  with  the  experimental  bound. 

As  just  described,  the  tensile  void  and  shear  band  failure  models,  when  used  in 
conjunction  with  the  slideline  erosion  algorithm  and  a  particular  mesh  for  the  jet 
particle  and  the  target  plate,  allowed  EPIC-2  to  produce  reasonable  agreement  with 
a  single  experiment  in  terms  of  final  hole  geometry,  the  rate  of  hole  formation,  and 
the  net  mass  lost  by  the  target  plate.  Three  problems  were  then  rerun  with  a  target 
mesh  having  double  the  fineness  it  had  in  all  previous  runs.  One  of  these  problems 
involved  no  failure  modeling,  one  involved  the  tensile  void  failure  model  only,  and 
the  third  involved  both  the  tensile  void  and  the  shear  band  failure  models.  For  the 
two  problems  involving  no  failure  modeling  or  only  tensile  failure  modeling,  no 
significant  changes  in  rho[e,  AM  and  r.95  were  produced  by  doubling  the  target  mesh 
fineness.  The  final  perforation  hole  remained  too  small  and  its  formation  completed 
prematurely  in  comparison  with  experiment.  The  problem  with  both  tensile  and 
shear  banding  failure  models  active  that  was  rerun  was  the  case  of  £perode  and  Aefau 
equal  to  0.50  and  0.00,  respectively.  This  is  the  case  that  with  the  original  mesh 
gave  best  overall  agreement  with  experiment  in  terms  of  rhoie,  AM,  t  95,  the  final 
hole  profile  and  shear  band  location  within  a  plate  cross-section.  Results  for  this 
problem  with  the  refined  mesh  showed  reductions  in  rhoU,  AM  and  t .95  from  values 
obtained  with  the  original  mesh  by  5%,  8.5%  and  13.7%,  respectively.  The  final 
hole  profile  became  flatter  (i.e.,  less  elevation  dependence  in  hole  radius)  than  with 
the  original  mesh.  The  location  of  shear  bands  also  became  less  dependent  on 
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elevation.  For  this  problem  with  both  failure  models  active,  the  CPU  time  required 
on  a  Cray  X-MP/48  demonstrated  a  factor  of  9.0  increase  over  the  time  required 
with  the  original  mesh.  In  light  of  this  increase,  the  solution  obtained  for  this 
problem  with  the  original  was  judged  to  be  sufficiently  accurate  to  be  useful. 

Possibilities  for  follow-up  work  will  now  be  considered.  In  the  present  study 
the  criteria  of  rhoie,  AM  and  t  95  agreement  with  experiment  have  not  uniquely 
determined  the  proper  set  of  input  values  for  the  three  parameters  e^ro*,  pfaii  and 
A zpfaU.  In  particular,  with  pfaU  fixed  at  -3.00  GPa  for  both  RHA  and  copper, 
agreement  was  obtained  for  a  range  of  Eperode  and  A Efait  combinations.  This 
consideration  motivated  the  recently  conducted  set  of  experiments  involving  various 
conditions  of  standoff  and  target  plate  thickness  (Raftenberg  1991,  1992).  The  same 
type  of  shaped  charge  warhead  used  in  the  two  experiments  cited  in  this  report  was 
also  employed  throughout  this  newer  set  of  experiments.  The  question  raised  is 
whether  all  experiments  in  this  combined  database  can  be  reasonably  well-matched 
in  terms  of  rhoie,  AM  and  t  95  by  a  single  set  of  input  values  for  parameters  e£r0(fc, 
Pfau  and  Atfait.  If  satisfactory  agreement  cannot  be  obtained  or  if  the  three 
parameters  remain  under-determined,  some  revision  of  the  material  models  used  in 
the  calculations  may  be  indicated.  These  would  include  the  models  for  flow  stress, 
tensile  failure  and  shear  banding  failure. 
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Table  1.  Parameter  Values  for  the  Leading  Jet 
Particle’s  Geometric  Model 


Parameter 

■MISRSTulll 

4.77 

12.75 

JMlf;  »n*llW 

Table  2.  Material  Parameter  Values  for  the  Metals 


Material 

Parameter 

OFHC 

CooDer 

RHA 

p0  (kg/m3) 

8945. 

7823. 

G  (GPa) 

46.31 

77.50 

Kx  (GPa) 

137.2 

163.9 

K 2  (GPa) 

175.1 

294.4 

K 3  (GPa) 

564.2 

500.0 

r 

1.96 

1.16 

A  (GPa) 

0.08963 

0.7922 

B  (GPa) 

0.2916 

0.5095 

N 

0.31 

0.26 

C 

0.025 

0.014 

M 

1.09 

1.03 

Tr  (K) 

294.3 

294.3 

1m  (K) _ 

1355.9 

1793.2 
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Table  3.  Results  from  Problem  Series  I 

(No  Failure  Modeling  in  RHA  or  Cu) 


Problem 

®  erode 

ffiole 

(mm) 

AM 

(8) 

*.95 

(Its) 

I-A 

0.25 

9.6 

27.1 

15.9 

I-B 

0.50 

10.3 

27.4 

16.4 

I-C 

0.75 

11.0 

27.7 

16.1 

I-D 

1.25 

12.1 

23.4 

15.6 

I-E 

2.00 

13.9 

20.3 

19.4 

Table  4.  Results  from  Problem  Series  II 

(Tensile  Failure  Modeling  in  RHA  &  Cu,  pfail  =  -  3.00  GPa; 
No  Shear  Band  Failure  Modeling) 


Problem 

E-erode 

rhole 

(mm) 

AM 

(8) 

*95 

(Its) 

II-A 

0.25 

9.7 

27.9 

10.5 

II-B 

0.50 

10.6 

28.6 

15.3 

II-C 

0.75 

10.9 

27.2 

22.7 

II-D 

1.25 

12.2 

24.9 

23.1 

II-E 

2.00 

15.0 

22.8 

89.3 
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Table  5.  Results  from  Problem  Series  III  Through  VIII 

(Tensile  Failure  Modeling  in  RHA  &  Cu,  pfaii  =  -  3.00  GPa; 
Shear  Band  Failure  Modeling  in  RHA,  A zfaU  Varied) 


Table  6.  Results  with  the  Refined  Mesh 


Problem 

£ erode 

A  £  fail 

Pfail 

(GPa) 

rhole 

(mm) 

AM 

(g) 

t.9S 

(Its) 

I-C(R) 

■ 

N.  A. 

N.  A. 

;';-V  '  ■ 

21.3 

II-C(R) 

N.  A. 

ETflEE 

ilKVSVH 

III-B(R) 

mmm 

0.00 

BaBlM 

BB1 

Table  7.  CPU  Times  on  a  Cray  X-MP/48  for  Six  Problems 


Problem 

£  erode 

^£-fail 

Pfail 

(GPa) 

CPU  Time 
(hrs) 

IQi 

■ 

N.  A. 

0.93 

UslfcM 

wmm 

N.  A. 

RES 

7.19 

TS— 

mausm 

N.  A. 

mgm 

0.93 

lllwal 

WSSM 

N.A. 

warn 

7.27 

ESI 

0.00 

mwm 

WESSM 

iS'i&ali  »Ta  »  mw4£l 
y  IIP  >11 

WMm 

0.00 

wmm 

mmsmi 
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RHA  PLATE 


X-RAY  TUBES  2000mm 
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43 


Figure  8.  Construction  of  Radial  Slices  to  be  Removed  From  Perforated  Plate. 
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Figure  10. 
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Problem  I-C  (No  Failure  in  RHA  &  Cu, 


lesh  Plot  at  4.02  us  After  Impact  for  Problem  I-C. 


Problem  I-C  (No  Failure  in  RHA  &  Cu,  tperode  -  0.7 


Problem  I-C, 


(WH) 
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NOTE:  r^ig  calculation  based  on  red  portion  of  boundary. 


Figure  20.  Mesh  Boundary  at  500.04  ps  and  Calculation  of  From  Problem  I-C. 
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Figure  21.  Computed  vs.  zperode  for  Problem  Series  I. 
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Rhole-avg  (mm) 


Problem  Series  1  (No  failure  in  RHA  or  Cu) 


Figure  23.  Mass  Lost  From  Target  Plate  vs.  Time  for  Problem  Series  I. 
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Figure  24.  Micrograph  of  a  Fracture  Surface  Formed  bv  Void 
Coalescence  in  the  RHA  Target  Plate  of  Rd.  10771 
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Optical  Micrograph,  800x,  2%  Nital  etch. 


Figure  26.  Micrograph  of  a  Crack  That  Extends  Along  a 
Shear  Band  in  the  RHA  Plate  of  Rd.  10771. 
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Figure  27.  Sketch  Showing  Locations  of  Prominent  Voids  and  Shear  Bands 
in  a  Radial  Slice  From  the  RHA  Plate  of  Rd.  10771. 
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Problem  Series  n  (Tensile  Failure  in  RHA  &  Cu,  PfaU  =  -3.00  GPa; 
No  Shear  Band  Failure) 
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Figure  35.  Computed  rhole  vs.  £%rode  for  Problem  Series  II. 
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Problem  Series  II  (Tensile  Failure  in  RHA  &  Cu,  pjaii  =  -3.00  GPa; 
No  Shear  Band  Failure) 


Figure  36.  Mass  Lost  From  Target  Plate  vs.  Time  for  Problem  Series  II. 
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Problem  Series  III  Through  VIII  (Tensile  Failure  in  RHA  &  Cu,  pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Varied) 
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Figure  38.  Computed  rhoie  vs.  £^0(jg  With  A Epfail  a  Parameter 
for  Problem  Series  III  Through  VIII. 
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(uiui)  3/°V  j 


Problem  Series  III  Through  VIII  (Tensile  Failure  in  RHA  &  Cu,  pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Ae^  Varied) 
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Problem 


have  failed  by  shear  banding. 


Problem  III-B  (Tensile  Failure  in  RHA  &  Cu,  pjait  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Ae?  ,  =  0.00;  tperode  =  0.50) 
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Figure  41.  Mesh  Plot  at  2.01  us  After  Impact  for  Problem  III-B. 


Problem  III-B  (Tensile  Failure  in  RHA  &  Cu,  pjail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  =  0.00;  zper0iu  =  0.50) 
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Figure  43.  Mesh  Plot  at  4.03  us  After  Impact  for  Problem  M-B 


Problem  III-B  (Tensile  Failure  in  RHA  &  Cu,  pjai!  =  -3.00  GPa 


Problem 


Problem  III-B  (Tensile  Failure  in  RHA  &  Cu,  Pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Ae^  =  0.00;  tper0<u  =  0.50) 


for  Problem 
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r  (fin) 

Figure  46.  Mesh  Plot  at  20.04  us  After  Impact  for  Problem  III-B. 


Problem  III-B  (Tensile  Failure  in  RHA  &  Cu,  pjaii  =  -3.00  GPa; 

■Shear  Band  Failure  in  RHA,  AeL,  =  0.00;  tperode  =  0.50) 
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Figure  47.  Mesh  Plot  at  50.04  us  After  Impact  for  Problem  III-B. 
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Problem  III-B  (Tensile  Failure  in  RHA  &  Cu,  Pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Ae^  =  0.00;  (Perode  =  0.50) 
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r  (hm) 

Figure  50.  Mesh  Plot  at  200.06  us  After  Impact  for  Problem  III-B. 


Problem  II1-B  (Tensile  Failure  in  RHA  &  Cu,  Pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  A£^  =  0.00;  tperode  =  0.50) 


Problem 


Problem  IV-C  (Tensile  Failure  in  RHA  &  Cu,  pfai}  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Ae^  =  0.10;  tperode  =  0.75) 


Problem  IV-C. 


Problem  V-D  (Tensile  Failure  in  RHA  &  Cu,  Pfmt  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Ae^ai/  =  0.25;  tper0tu  =  1-25) 
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Figure  56.  Mesh  Plot  at  500.06  us  After  Impact  for  Problem  V-D. 


Problem  VI-E  (Tensile  Failure  in  RHA  &  Cu,  pj-aii  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  A eL/  =  0.50;  tperode  =  2.00) 
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Figure  57.  Mesh  Plot  at  500.04  us  After  Impact  for  Problem  VI-E. 
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Shear  Band  Failure  in  RHA,  A ^  =  0.50) 
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Problem  Series  VII  (Tensile  Failure  in  RHA  &  Cu,  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  A ^  =  0.75) 
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Problem  Series  VIII  (Tensile  Failure  in  RHA  &  Cu,  pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Ae^  =  1.00) 


Figure  63.  Mass  Lost  From  Target  Plate  vs.  Time  for  Problem  Series  VIII. 
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(3)  JVV 


Problem  Series  m  Through  VIII  (Tensile  Failure  in  RHA  &  Cu,  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  £^ail  Varied) 
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Figure  65.  Computed  AM  vs.  A e&a  With  zpprnde  a  Parameter 
for  Problem  Series  III  Through  VIII. 
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Figure  66.  The  Initial  Mesh  Used  in  Section  9. 


Problem  I-C(R)  (No  Failure  in  RHA  or  CU;  tperode  =  0.75; 
Refined  Target  Mesh) 
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Figure  69.  Mesh  Plot  at  500.02  us  After  Impact  for  Problem  I-C(R). 
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Problem  II-C(R)  (Tensile  Failure  in  RHA  &  Cu,  pjail  =  -3.00  GPa; 
No  Shear  Band  Failure;  =  0.75; 

Refined  Target  Mesh) 
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Figure  71.  Mesh  Plot  at  1,00  us  After  Impact  for  Problem  II-C(R). 


Problem  II-C(R)  (Tensile  Failure  in  RHA  &  Cu,  pjail  =  -3.00  GPa; 
No  Shear  Band  Failure;  e^rode  =  0.75; 

Refined  Target  Mesh) 
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Figure  72.  Mesh  Plot  at  3.01  us  After  Impact  for  Problem  H-C(R). 


Problem  H-C(R)  (Tensile  Failure  in  RHA  &  Cu,  Pfaj  =  -3.00  GPa; 
No  Shear  Band  Failure;  £%rodc  =  0.75; 

Refined  Target  Mesh) 
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Figure  73.  Mesh  Plot  at  500.00  us  After  Impact  for  Problem  D-C(R). 
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Problem  III-B(R)  (Tensile  Failure  in  RHA  &  Cu,  pfat!  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  A£jaU  =  0.00;  efroc^  =  0.50; 
Refined  Target  Mesh) 
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Figure  76.  Mesh  Plot  at  2.01  us  After  Impact  for  Problem  III-B(R). 


Problem  III-B(R)  (Tensile  Failure  in  RHA  &  Cu,  pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Aj$ail  =  0.00;  zprode  =  0.50; 
Refined  Target  Mesh) 
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Figure  77.  Mesh  Plot  at  3.01  us  After  Impact  for  Problem  III-B(R). 


Problem  III-B(R)  (Tensile  Failure  in  RHA  &  Cu,  pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  te?faiS  =  0.00;  tperode  =  0.50; 
Refined  Target  Mesh) 
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Figure  78.  Mesh  Plot  at  4.00  us  After  Impact  for  Problem  DI-B(R). 


Problem  III-B(R)  (Tensile  Failure  in  RHA  &  Cu,  pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  Ae^  =  0.00;  tperode  =  0.50; 
Refined  Target  Mesh) 
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Figure  80.  Mesh  Plot  at  10.02  us  After  Impact  for  Problem  III-B(R). 
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Shear  Band  Failure  in  RHA,  Ae^ai/  =  0.00;  tperodt  =  0.50; 
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Figure  82.  Mesh  Plot  at  50.03  us  After  Impact  for  Problem  III-B(R). 
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Figure  83.  Mesh  Plot  at  100.01  us  After  Impact  for  Problem  III-B(R'). 
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Figure  85.  Mesh  Plot  at  500.01  us  After  Impact  for  Problem  ni-B(R). 
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Piobicin  in-B(R)  (Tensile  Failure  in  RHA  &  Cu,  Pfail  =  -3.00  GPa; 

Shear  Band  Failure  in  RHA,  A^ail  =  0.00;  tperode  =  0.50; 
Refined  Target  Mesh) 


2  r  cycle  =  14907  time  =  500.01  ps 


NOTE:  rfofe  calculation  based  on  red  portion  of  boundary. 

Figure  86.  Mesh  Boundary  at  500.01  ps  and  Calculation  of  rtote  for  Problem  m-B(R). 
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Figure  87. 
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LIST  OF  SYMBOLS 


A,  B,  C,  M,  N 
Dp 

G 

Ku  K2,  K 3 
Lc,  Ls,  Lb,  Rs,  Rb 


Lp 


T 

Tm 

Tr 

T* 

VP 

e 

f{ep) 

P 

Pfail 


rhole 


S rri  $ZZ’  ^68 


parameters  in  the  Johnson-Cook  strength  model 

ij  component  of  the  plastic  part  of  the  rate  of  deformation 
tensor 

elastic  shear  modulus 

parameters  in  the  Mie-Gruneison  equation  of  state 

parameters  in  the  geometric  model  for  the  leading  jet  particle 

length  of  the  leading  jet  particle 

temperature 

melting  temperature 

room  temperature 

homologous  temperature 

volume  of  the  leading  jet  particle 

internal  energy  per  underformed  volume 

stress  reduction  function  used  in  the  shear  banding  failure 
model 

pressure 

pressure  at  which  the  onset  of  tensile  failure  occurs 
radial  coordinate 

computational  hole  radius  at  500  its  after  initial  impact, 
averaged  over  the  target  plate’s  thickness 

experimental  final  hole  radius  averaged  over  the  target  plate’s 
thickness 

ij  component  of  the  deviatoric  Cauchy  stress  tensor  following 
reduction 

radial,  axial,  circumferential  normal  components  of  the 
deviatoric  Cauchy  stress  tensor  following  reduction 
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j 


sn  shear  component  of  the  (deviatoric)  Cauchy  stress  tensor 

following  reduction 

sjj  ij  component  of  the  deviatoric  Cauchy  stress  tensor  prior  to 

reduction 

s *r,  szz,  sqq  radial,  axial,  circumferential  normal  components  of  the 

deviatoric  Cauchy  stress  tensor  prior  to  reduction 

s*z  shear  component  of  the  (deviatoric)  Cauchy  stress  tensor  prior 

to  reduction 

t  time 

t. 95  computational  time,  after  initial  impact,  at  which  the  target 

plate’s  cumulative  mass  loss  equals  0.95AM 

t. 9T  experimental  time,  after  initial  impact,  at  which  the  target 

plate’s  cumulative  mass  loss  equals  0.95A Mexp- 

z  axial  coordinate 

zp  axial  coordinate  of  the  centroid  of  the  leading  jet  particle 

AM  computational  accumulated  mass  lost  by  the  target  plate  at  500 

Its  after  initial  impact 

A Mexp  experimental  total  mass  lost  by  the  target  plate 

A Epfau  the  minimum  increment  above  Eponset  of  equivalent  plastic  strain 

corresponding  to  which  all  components  of  deviatoric  Cauchy 
stress  within  an  element  are  set  to  zero 

T  Gruneison  coefficient 

X  expression  involving  parameters  of  the  geometric  model  for 

the  leading  jet  particle 

tp  equivalent  plastic  strain 

£Prode  equivalent  plastic  strain  at  which  a  slideline’s  master  surface 

element  is  eroded 

EPonse,  equivalent  plastic  strain  within  an  element  at  the  time  of  the 

onset  of  its  failure  by  shear  banding 
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time  rate  of  change  of  equivalent  plastic  strain 
fractional  change  in  density  from  the  undeformed  state 
density 

undeformed  density 
equivalent  Cauchy  stress 
flow  stress 

radial,  axial,  circumferential  normal  components  of  the  Cauchy 


stress  tensor 


Intentionally  left  blank. 


130 


No.  of 
Copies 

2 

1 

1 

2 

2 

1 


(Uncials,  only)  } 

1 


Organization 


No.  of 

Copies  Organization 


Administrator 

Defense  Technical  Info  Center 
ATTN:  DliC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCAM 
5001  Eisenhower  Ave. 

Alexandria,  VA  22333-0001 

Commander 

U.S.  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
2800  Powder  Mill  Rd. 

Adelphi,  MD  20783-1145 

Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-IMI-I 
Pica  tinny  Arsenal,  NJ  07806-5000 

Commander 

U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 
ATTN:  SMCAR-TDC 
Pica  tinny  Arsenal,  NJ  07806-5000 

Director 

Benet  Weapons  Laboratory 
U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  1218<M050 

Commander 

U.S.  Army  Rock  Island  Arsenal 
ATTN:  SMCRI-TL/Technical  Library 
Rock  Island,  IL  61299-5000 

Director 

U.S.  Army  Aviation  Research 
and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 

Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOQ 
Redstone  Arsenal,  AL  35898-5010 


1  Commander 

U.S.  Army  Tank-Automotive  Command 
ATTN:  ASQNC-TAC-DIT  (Technical 
Information  Center) 

Warren,  MI  48397-5000 

1  Director 

U.S.  Army  TRADOC  Analysis  Command 
ATTN:  ATRC-WSR 

White  Sands  Missile  Range,  NM  88002-5502 

1  Commandant 

U.S.  Army  Field  Artillery  School 
ATTN:  ATSF-CSI 
Ft  Sill,  OK  73503-5000 


(cum.  only)  i  Commandant 

U.S.  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 
Fort  Benning,  GA  31905-5660 

(UbcUm.  oily)  i  Commandant 

U.S.  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 

1  WL/MNOI 

Eglin  AFB,  FL  32542-5000 

Aberdeen  Proving  Ground 

2  Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 


1  Cdr,  USATECOM 
ATTN:  AMSTE-TC 

3  Cdr,  CRDEC,  AMCCOM 
ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSI 

1  Dir,  VLAMO 

ATTN:  AMSLC-VL-D 

10  Dir,  USABRL 

ATTN:  SLCBR-DD-T 


131 


No.  of 

Copies  Organization 


No.  of 

Copies  Organization 


2  Director 

Defense  Advanced  Research  Projects 
Agency 

ATTN:  J.  Richardson 

Technical  Information 
1400  Wilson  Blvd. 

Arlington,  VA  22209-2308 

3  U.S.  Army  Research  Office 
ATTN:  J.  Chandra 

K.  Iyer 

Technical  Library 
P.O.  Box  12211 
4300  Miami  Blvd. 

Research  Triangle  Park,  NC  27709 

1  U.S.  Army  MICOM 

ATTN:  AMSMI-RD-ST-WF, 

Michael  Schexnayder 
Redstone  Arsenal,  AL  35898-5250 

5  U.S.  Army  Materials  Technology  Laboratory 
ATTN:  SLCMT-MRT, 

Tony  Chou 
Charles  S.  White 
A.M.  Rajendran 
J.  Beatty 

Technical  Library 
Watertown,  MA  02172-0001 

1  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-AA,  J.  Pearson 
Picatinny  Arsenal,  Ni  07806-5000 

2  Commander 

Naval  Surface  Warfare  Center 
ATIN:  W.  Holt 

Technical  Library 
Dahlgren,  VA  22448-5000 

4  Commander 

Naval  Surface  Warfare  Center 
ATTN:  F.  J.  Zerilli 

Code  U12,  H.  Chen 
Code  R14,  H.  Mair 
Code  X211,  Library 
Silver  Spring,  MD  20902-5000 


1  U.S.  Naval  Academy 
Dept  of  Mathematics 
ATTN:  R.  Malek-Madani 
Annapolis,  MD  21402 

2  Naval  Weapons  Center 
ATTN:  S.  A.  Finnegan 

Technical  Library 
China  Lake,  CA  93555 

3  Air  Force  Armament  Laboratory 
ATTN:  J.  Foster 

W.  Cook 
Technical  Library 
Eglin  AFB,  FL  32542-5438 

1  Air  Force  Wright  Aeronautical  Laboratories 

Air  Force  Systems  Command  Materials 
Laboratory 
ATTN:  T.  Nicholas 
Wright-Patlerson  AFB,  OH  45433 

9  Director 

Los  Alamos  National  Laboratory 
ATTN:  G.  E.  Con,  F663 

Thomas  F.  Adams,  F663 
Robert  Karpp,  MS  J960 
John  Bolstad,  MS  G787 
K.  S.  Holian,  MS  B295 
J.  Chapyak 
P.  Follansbee 
G.  T.  Gray 
Technical  Library 
P.O.  Box  1663 
Los  Alamos,  NM  87545 

1  Sandia  National  Laboratories 

ATTN:  D.  Bammann 
Livermore,  CA  94550 

1  University  of  Missouri-Rolla 

Dept  of  Mechanical  and 
Aerospace  Engineering 
ATTN:  R.  Batra 
Rolla,  MO  65401-0249 


132 


No.  of 

Copies  Organization 


No.  of 

Copies  Organization 


9  Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  Barry  R.  Bowman,  L-122 
John  E.  Reaugh,  L-290 
Michael  J.  Murphy 
W.  H.  Gourdin 
R.  Whirley 
G.  Goudreau 
D.  J.  Steinberg,  L-35 
Robert  Tipton,  L-35 
Technical  Library 
P.O.  Box  808 
Livermore,  CA  94550 

7  Sandia  National  Laboratories 

ATTN:  L.  Davison 
W.  Herman 
M.  Forrestal 
J.  M.  McGlaun 
P.  Yarrington 
V.  Luk 
D.  Grady 
P.O.  Box  5800 

Albuquerque,  NM  87185-5800 

1  Denver  Research  Institute  University 

of  Denver 
ATTN:  R.  Recht 
P.O.  Box  10127 
Denver,  CO  80210 

1  Massachusetts  Institute  of  Technology 

Dept,  of  Mechanical  Engineering 
ATTN:  L.  Anand 
Cambridge,  MA  02139 

1  California  Research  and  Techology,  Inc. 

ATTN:  Marie  Majerus 
P.O.  Box  2229 
Princeton,  NJ  08543-2229 

1  Rennselaer  Polytechnic  Institute 

Dept,  of  Computer  Science 
ATTN:  J.  Flaherty 
Troy,  NY  12181 


3  Brown  University 

Division  of  Engineering 
ATTN:  R.  Clifton 

J.  Duffy 

A.  Needle  man 
Providence,  RI  02912 

1  Cornell  University 

Dept,  of  Theoretical  and  Applied  Mechanics 

ATIN:  J.  Jenkins 
Ithaca,  NY  14850 

1  North  Carolina  State  University 

Dept,  of  Civil  Engineering 
ATTN:  Y.  Horie 
Raleigh,  NC  27607 

1  North  Carolina  State  University 

Dept  of  Mechanical  Engineering 

ATTN:  M.  Zikry 

Raleigh,  NC  27607 

2  The  Johns  Hopkins  University 
Dept  of  Mechanical  Engineering 
Latrobe  Hall 

ATTN:  A.  Douglas 

K.  T.  Ramesh 
34th  and  Charles  Streets 
Baltimore,  MD  21218 

3  University  of  California 

at  San  Diego 

DepL  of  Applied  Mechanics  and 
Engineering  Sciences 
ATTN:  S.  Nemat-Nasser 
M.  Meyers 
D.  J.  Benson 
La  Jolla,  CA  92093 

1  Northwestern  University 

Dept  of  Applied  Mathematics 
ATTN:  T.  Belytschko 
Evanston,  IL  60201 

1  University  of  Illinois 

Dept  of  Theoretical  and 
Applied  Mechanics 
ATTN:  T.  Shawki 
Urbana,  IL  61801 


133 


No.  of 

Copies  Organization 

1  University  of  Texas 

Dept,  of  Engineering  Mechanics 
ATTN:  J.  Oden 
Austin,  TX  78712 

1  University  of  Maryland 

Dept  of  Mechanical  Engineering 
ATTN:  R.  Armstrong 
College  Park,  MD  20742 

1  Washington  State  University 

Dept,  of  Physics 
ATTN:  Y.  Gupta 
Pullman,  WA  99163 

1  Washington  State  University 

Dept  of  Mechanical  and  Materials 
Engineering 
ATTN:  H.  Zbib 
Pullman,  WA  99164-2920 

1  Arizona  State  University 

Mechanical  and  Aerospace  Engineering 
M/S  6106 

ATTN:  D.  Krajcinovic 
Tempe.AZ  85287-6106 

1  Michigan  Technological  University 

Dept  of  Mechanical  Engineering 
ATTN:  William  W.  Predebon 
1400  Townsend  Dr. 

Houghton,  MI  49931-1295 

1  University  of  Alabama 

ATTN:  S.  Jones 
210  Hardaway  Hall 
Box  870278 

Tuscaloosa,  AL  35487-0278 

1  University  of  Texas  at  El  Paso 

Dept  of  Metallurgical  and  Materials 
Engineering 
ATTN:  L.  E.  Murr 
El  Paso,  TX  79968 


No.  of 

Copies  Organization 

1  State  University  of  New  York  at  Stony  Brook 
Applied  Mathematics  and  Statistics  Department 
ATTN:  James  Glimm 
Stony  Brook,  NY  11794 

1  Aerojet  Precision  Weapons 
Dept.  5131/T-W 

ATTN:  Joseph  Carleone 
1 100  Hollyvale 
Azusa,  CA  91702 

2  Dyna  East  Corporation 
ATTN:  P.  C.  Chou 

W.  Flis 
3201  Arch  St 
Philadelphia,  PA  19104 

2  Southwest  Research  Institute 
ATTN:  C.  Anderson 
S.  Bodner 
P.O.  Drawer  28255 
San  Antonio,  TX  78228-0255 

1  Battelle 
Edgewood  Operations 
ATTN:  R.  Jameson 
2113  Emmorton  Park  Rd. 

Edgewood,  MD  21040 

2  Alliant  Techsystems,  Inc. 

ATTN:  Gordon  R.  Johnson 

Tim  Holmquist 
MN  48-2700 
7225  Northland  Dr. 

Brooklyn  Park,  MN  55428 

1  S-Cubed 

ATTN:  R.  Sedgwick 
P.O.  Box  1620 
La  Jolla,  CA  92038-1620 

1  Orlando  Technology,  Inc. 

ATTN:  Dan  Matuslca 
P.O.  Box  855 
Shalimar,  FL  32579 
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1  Livermore  Software  Technology  Corp. 
ATTN:  John  O.  Hallquist 

2876  Waverly  Way 
Livermore,  CA  94550 

2  SRI  International 
ATTN:  D.  Curran 

H.  Giovanola 
333  Ravens  wood  Ave. 

Menlo  Park,  CA  94025 
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1  Defense  Research  Suffield 

ATTN:  Chris  Weickert 
Ralston,  Alberta,  JOJ  2NO  Ralston 
Canada 

1  Ernst  Mach  Institute 

ATTN:  A.  J.  Stilp 
Eckerstrasse  4 
D-7800  Frieburg  i.  Br. 

Germany 

1  Royal  Armament  R&D  Establishment 

ATTN:  Ian  Cullis 
Fort  Halstead 

Sevenoaks,  Kent  TN14  7BJ 
England 

1  Centre  d’Etudes  de  Gramat 

ATTN:  SOLVE  Gerard 
46500  Gramat 
France 

1  Materials  Research  Laboratories 

ATTN:  R.  L.  Woodward 
P.O.  Box  50 
Ascot  Vale 
Victoria  3032 
Australia 

1  Israel  Military  Industries 

ATTN:  E.  Hirsch 
Ramat  Hasharon  47100 
Israel 
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